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Abstract 

We develop an original approach for the quantitative characterisation of the en- 
tanglement properties of, possibly mixed, bi- and multipartite quantum states of 
arbitrary finite dimension. Particular emphasis is given to the derivation of reliable 
estimates which allow for an efficient evaluation of a specific entanglement measure, 
concurrence, for further implementation in the monitoring of the time evolution of 
multipartite entanglement under incoherent environment coupling. The flexibility 
of the technical machinery established here is illustrated by its implementation for 
different, realistic experimental scenarios. 
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1 Introduction 



Entanglement is one of the central issues of debate in quantum theory since 
the beginning of the last century and certainly a key idea when it comes to dis- 
tinguish classical and quantum concepts. Moreover, besides this fundamental 
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aspect, the interest in entangled states has been recently renewed because their 
properties lie at the heart of many potential applications. Be it in quantum 
computation [1,2], teleportation [3] or quantum cryptography [4], entangle- 
ment is viewed as an important resource and, as such, must be quantified. In 
addition, great experimental progresses in the production, manipulation and 
detection of entangled states [5,6,7,8,9,10,11,12,13,14,15,16,17] require such a 
quantification to be versatile enough to deal with the states encountered in 
actual experiments, which are in general mixed and typically involve several 
particles. 

The first attempt to discern the non-local correlations of measurement re- 
sults induced by entanglement was formulated with Bell's inequalities [18,19], 
which underwent a first experimental check [20] in the sixties. Bell's inequal- 
ities are capable of discriminating correlations due to entanglement against 
those described by local hidden variable models [21]. Later, also entanglement 
criteria that use special three-partite states, without involving inequalities, 
were found [22] and tested experimentally [7] . Albeit able to reveal the entan- 
gled nature of some quantum states, the above criteria cannot (and do not 
intend) to quantify the amount of entanglement carried by a given state. 

Only recently has the problem of finding a quantity that measures quan- 
tum correlations been studied more intensively [23,24,25]. Virtually the entire 
state-of-the-art theory of entangled quantum states is based on so-called en- 
tanglement measures, scalar quantities that quantify quantum correlations, 
and distinguish them from classical ones. For bipartite pure states such mea- 
sures exist and are straightforwardly computable. However, if one aspires to 
describe realistic states observed in experiments, it is imperative to allow for 
a proper quantification also of mixed states entanglement, since there is no 
system that could be decoupled perfectly from environmental influences, and 
mixing is thus unavoidable. 

Although several measures for mixed states have been proposed, no simple 
criterion of discriminating classical from quantum correlations is known so 
far. All proposed measures that unambiguously fulfill this task involve some 
- generally high dimensional - optimisation procedure, and hardly allow for 
an explicit evaluation in concrete cases. 1 Only for states of smallest possible 
dimensions, i.e., for bipartite two-level systems, do a few measures exist that 
can be evaluated algebraically [27,28,29], though no such measures are known 
for arbitrary states of higher dimensional systems. 

Hence, there is some mismatch between the more formal part of entangle- 
ment theory, which seeks for the general characterisation of arbitrary quantum 
states, and experimental progresses of the last years, which are leading to the 

1 Negativity [26] can be evaluated algebraically, though does not detect all entan- 
gled states. 
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production of (possibly specific classes of) ever more complex entangled states 
of multipartite quantum systems. For a theoretical analysis of the latter, we 
need algebraic as well as numerical tools to describe static as well as dynamical 
properties of multipartite quantum states in a quantitative manner - what is 
a nontrivial target, due to the rapidly increasing complexity of the problem 
at hand, as particle number and effective Hilbert space dimensions are in- 
creased. Studying the dynamics of entanglement requires efficiently evaluable 
quantifiers thereof as an undispensible ingredient, ideally for arbitrary system 
sizes. 

In the present review, we attempt to improve on this mismatch, by responding 
precisely to this latter requirement: Starting from the formal algebraic descrip- 
tion of a suitable entanglement measure - concurrence -, we derive a hierarchy 
of bounds [30] and approximations [31] thereof, which imply progressively re- 
duced computational efforts for its actual evaluation for bi- or multipartite, 
mixed quantum states in arbitrary finite dimensions. After numerical tests of 
the tightness of these various estimates, we implement this novel toolbox, to 
monitor entanglement dynamics under experimentally realistic conditions [32]. 

We start out with a brief recollection of the most important concepts of en- 
tanglement theory. 

1.1 Entangled states 

1.1.1 Bipartite entanglement 

A bipartite system is a quantum system that is composed of two physically 
distinct subsystems. It is associated with a Hilbert space Ti that is given by 
the tensor product Tii ® Ti.2 of the predefined factor spaces, each of which 
describes one subsystem. 

For pure states one distinguishes two different kinds of states. A state is 
called a product state or separable, if it can be written as a tensor product of 
subsystem states, i.e. if there are states \<p) and \<p) of the subsystems such 
that 



Such a state describes a situation analogous to a classical one insofar as the 
system state contains exactly the information that is contained in the subsys- 
tem states. A state reduction [33] caused by a measurement performed on one 
subsystem has no influence on the state of the other subsystem. This means 
that measurement results on the different subsystems are uncorrelated (or in- 
dependent). In contrast to this, they are correlated for entangled states, i.e. 
states that cannot be written as a product of subsystem states as in Eq. (1). 
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Here, a local measurement causes a state reduction of the entire, i.e. bipar- 
tite, system state, and therefore changes the probabilities for potential future 
measurements on either subsystem. 

For mixed states the situation is more complicated. Product and separable 
states are not synonymous anymore. Whereas the former can be expressed as 
a tensor product of subsystem states, i.e. there are states Q\ and Q2 describing 
the single subsystems, such that 

Q = Qi ® 02 , (2) 

a convex sum of such product states is needed to represent a general separable 
state 

q = E Pi q\ ® A , ( 3 ) 

i 

where convexity implies positive coefficients pi that sum up to unity, Y,i Pi — 1- 
Such a state refers to a situation where correlations between different sub- 
systems are due to incomplete knowledge about the system state. They are 
characterised completely by the classical probabilities Pi. 

Quantum correlations, i.e. , entanglement, need now be distinguished from 
classical correlations - a problem which we will focus on throughout the largest 
part of this paper. Formally, and in a rather non-constructive way, an entan- 
gled mixed state is defined through the non-existence [34] of a convex decom- 
position alike Eq. (3). 

Q ^ Pi e\ ® q\ , (4) 

i 

what is the mixed state generalization of the negation of (1). The correlations 
contained in such states cannot be characterised completely by a set of classical 
probabilities. Thus, entangled states bear correlations that do not exist in any 
classical system. 

1.1.2 Multipartite entanglement 

The definition of entangled states can straightforwardly be generalised to mul- 
tipartite systems, i. e. systems that decompose into more than two subsystems. 
A p-partite system is described by a Hilbert space 7i that decomposes into 
the tensor product of p factor spaces Tii <E> • • • <8> 7~L p . A pure state is separable 
if it can be written as a product of p states, each of which describes one of the 
subsystems - any state that is not separable is entangled. A mixed state is sep- 
arable if it can be written as a convex sum of product states, i. e. of products 
of states each acting on a single subsystem. If it cannot, it is entangled. 

For multipartite systems it is furthermore meaningful to distinguish between 
different degrees of separability: a pure state \^/) is called k-separable if it can 
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be written as a product of k states \<f>i), each of which is an element of one 
of the factor spaces, or of the product of some of those. Thus a ^-separable 
state with k < p can contain entanglement between some of the subsystems, 
whilst there also are subsystems that are completely uncorrected. In this 
terminology, p-separability is equivalent to complete separability. 

1.2 Separability criteria 

1.2.1 Pure state separability - The Schmidt decomposition 

Pure bipartite states can be classified with the help of their Schmidt decompo- 
sition. Each bipartite pure state can be expressed in some product basis, 

l*> = E ba Wi) ® l&> • (5) 

ij 

The local bases {| and can be chosen arbitrarily. However, refer- 

ring to a given state, there is always one distinguished basis. It can be con- 
structed with the following representations of the identity operators: l x = 
Yl,iU '\ipi)(ipi\U acting on the first factor space Tii, and, analogously I2 = 
J2i V^\(f>i){<f>i\V, acting on the second factor space H.2- U and V are some ar- 
bitrary, local unitary transformations on r hi 1 and Ti 2l respectively. Inserting 
these identities in Eq. (5), the state |\&) can be expressed as 

|¥) = £[uHy u^ t )®v^ 3 ) , (6) 

where the unitary matrices u and v are defined as 

Uij = ((pi\U\ipj) and Vij = {<f>j\V\<t>i) ■ 2 (7) 

Now one can use the fact that every complex matrix b can be diagonalised by 
two unitary transformations u and v, such that ubv, with real and non-negative 
diagonal elements Si, provides the singular value decomposition of b [35]. 
Hence, any pure state can be represented in terms of its Schmidt coefficients 
Aj = Sf, and of the associated Schmidt basis ® |£j) 2 = U^\<fi) <8> V^\<pi): 

l*> = E l&>i ® l&> 2 , (8) 

i 

with the sum limited by d, the dimension of the smaller subsystem. Given 
that the Schmidt basis comprises - by construction - only separable states, all 
information about the entanglement of |^) is now contained in the Schmidt 
coefficients Aj. The characterisation of all correlations of a given pure state is 

2 Note that the order of indices in the definition of v is different from that of u. 
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therefore tantamount to the knowledge of all Schmidt coefficients. The nor- 
malisation condition = 1 implies that there are d — 1 independent 
coefficients. 

The Schmidt coefficients can be easily computed with the help of one of the 
reduced density matrices Q\ = Tr2 1 (^\, Q2 = Tri|^)(^|. Assume, without 
loss of generality, d = dim(7ii) < dim(7Y 2 )- Using Eq. (8), one easily verifies 
that the spectrum of g\ is just given by the Schmidt coefficients - the spectrum 
of Q2 is given by the Schmidt coefficients and dim(7i2) — dim(7ii) vanishing 
eigenvalues. 

The Schmidt coefficients also allow to distinguish separable from entangled 
states - a separable state is characterised by a vector of Schmidt coefficients 
with only one non- vanishing entry: A = \ s = [1, 0, . . . , 0], whereas the Schmidt 
vector of an entangled state has at least two non-vanishing components. A 
state is called maximally entangled, if its Schmidt vector reads A = A m = 
[l/d, . . . , 1/d]. In Section 1.4 we will discuss in which respect this terminology 
is legitimate. 

It follows that the concept of Schmidt coefficients allows to relate the de- 
gree of entanglement of pure bipartite states to the degree of mixing of the 
corresponding reduced density matrices - a pure reduced density matrix cor- 
responds to a separable state, whereas a maximally entangled state leads to a 
maximally mixed reduced density matrix. 

1.2.2 Mixed state separability 

Whereas we have just seen that the separability of pure bipartite states can 
be easily checked, it turns out to be much more difficult to decide whether a 
given mixed state bears quantum entanglement. The above definition, Eq. (4), 
for entangled mixed states is not constructive, and generically it is not clear 
whether there is a set of product states such that g can be represented as a 
convex sum of its elements. 

The standard approach to decide on the separability of a given mixed state 
relies on positive maps. A map A : B(H) — > B(H), where 13(H) is the space of 
bounded linear operators on H, is called positive if it maps positive operators 
on positive ones, i.e. 

A(g) > 0, for all g > , (9) 

where positivity of an operator g is just a short hand notation stating that 
g is positive semi-definite, i.e. it has only non-negative eigenvalues. A crucial 
property of positive maps is that a trivial extension A <g> 1 is not necessarily 
positive [36]. Consider a positive map A : B(Hi) 1— > B(Hi): if the trivial 
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extension A <g> 1, with 1 the identity map on B(H,2), is not positive, it can be 
used to conclude on the separability of a mixed state g, acting on Tii ® Ti^- 
Since the extended map A ® 1 is not positive, there are some states rj such 
that (A ® l)(r?) ^ 0. However, if one assumes the considered state g to be 
separable, its convex decomposition into product states (3) implies 

(A®l)te)=£ Pi A(Q\)®Qi. (10) 

i 

Obviously, any expectation value of this quantity is non-negative, such that 
(A <g> l^j(g) > 0. Equivalently, a state g is entangled if (A (g) l)(£>) ^ 0. 

However, the inverse statement does not hold in general. The mere fact that 
g remains positive under the extended map does not necessarily imply that g 
is separable. Only if 

(A <g> l) (g) > for all positive maps A , (11) 

one may conclude that g is separable [37]. Note that for the complementary 
implication alone one only needs to find one positive map A with (A®lj(g) ^ 
0. This statement does not allow to derive a sufficient separability criterion 
for the very general case, since the classification of positive maps is still an 
unsolved problem. A large 3 class of entangled states is detected by the special 
choice of the transposition T =: A [38] - that indeed is a positive map. The 
partial transpose gP 1 = (T (g) lj(g) of a state g is deduced as the relevant 
auxiliary quantity: if g pt has at least one negative eigenvalue, 

f t , (12) 

the state g is entangled. 

However, if (f 1 is positive, one can infer separability of g only for low- dimensional, 
namely 2x2 and 2x3 systems. For these, the positive partial transpose (ppt) 
or standard criterion unambiguously distinguishes separable and entangled 
states [37]. However, in higher dimensions there exist entangled states [39,40] 
that are not detected by the ppt criterion. 

1.3 How to quantify entanglement? 

Since the definition of entangled states given in Eq. (4) is not constructive, it 
turned out difficult to decide whether a given state is separable - and the gen- 
eral solution of this problem is still unknown. Moreover, the non-constructive 

3 'large' in the topological sense: the set of states with negative partial transpose 
contains an open set. 
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definition also complicates finding a quantitative description of entanglement, 
rather than the purely qualitative one. How can you measure something, if 
you don't even know what it is? 

The basic idea for a quantitative treatment is to classify all kinds of operations 
that in principle can be applied to quantum systems and that can create or 
increase only classical correlations, but none of quantum nature. Any quantity 
that is supposed to quantify entanglement needs to be monotonously decreas- 
ing under such operations [25,41]. 

In our subsequent discussion, we will not distinguish between operations de- 
scribing the time evolution of a real system, and those which serve just as a 
mathematical tool. In the latter case one can always have in mind a Gedanken- 
experiment where the considered operation is implemented. For the following 
considerations it is not crucial whether one has the technical prerequisites and 
experimental skills to perform a considered operation - but rather that the 
operation is in principle allowed by the laws of quantum mechanics. There- 
fore, we do not consider technical problems - as long as we do not refer to real 
experiments. 

A map £ : B(Hi) — > B(Hf) describing the evolution of a quantum system has 
to be linear, 

£(\g + X'g') — A £(g) + X' £{q') , (13) 

due to the underlying linear Schrodinger equation. Moreover, in order to en- 
sure positivity of g, any map £ has to be positive. However, this requirement 
is not strong enough to ensure positivity of g in all cases. Since one can always 
consider a system as a subsystem of a larger one, one has to allow for exten- 
sions £ <E> 1 of £. The extended map acts on the entire system in such a way 
that the original map affects the considered subsystem, whereas the identity 
map acts on the residual system degrees of freedom. As already mentioned 
in Section 1.2.2, such a trivial extension does not necessarily preserve posi- 
tivity. In order to guarantee positivity of the entire system state, one has to 
require that the described extension be a positive map for identity maps 1 in 
any dimension, i. e. that £ is completely positive. Consequently, any evolution 
consistent with the general rules of quantum mechanics can be described by 
a linear, completely positive map, called quantum operation. 

A unitary evolution is just a special case of such a quantum operation - gen- 
eral quantum operations can also describe non-unitary evolutions, e.g. due to 
environment coupling or measurements. Any such quantum operation can be 
composed from elementary operations [36,42,43]: 

- unitary transformations, £\{q) = UqU^ ; 

- addition of an auxiliary system, £2(0) = £?®cx, where g is the original system 
and a is the auxiliary one; 
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- partial traces, 83(g) = Tr p g ; 

- projective measurements, 84(g) = PkgPk/ r £r:(P k g) , with P| = P k , 

which allows for a physical interpretation thereof. For a formal mathematical 
treatment it is useful to note that any quantum operation can be expressed 
as an operator sum [44,45] 

8(g) = £ E lQ E\ , with ]T e\e, = 1 , (14) 

i i 

with suitably defined linear operators Ei. 

The reduced dynamics of a system initially prepared in the state g(0), coupled 
to an environment with initial state |eo) can be interpreted in terms of quan- 
tum operations. If we allow for an interaction between system and environ- 
ment, we will have a unitary evolution U(t) in both system and environment. 
The system state after time t is obtained by evolving the system-bath state 
over t, followed by a trace over the environment: 

g(t) = Tr e nv( U(t) | e )(e | ® g(0) U ] (t) ) . (15) 

Expressing the trace over the environmental degrees of freedom by a sum over 
an orthonormal basis {|x?)}, one immediately obtains the above operator sum 
representation 

g(t) = ]T Ei(t) g(0) Ej(t) , with E % = ( Xi \U(t)\e Q ) , (16) 

i 

where the operators Ei(t) satisfy the resolution of the identity required in 
Eq. (14). 

For our purposes, it will be useful to distinguish the following types of quantum 
operations: 

- local operations, 

- global operations, 

- local operations and classical communication (LOCC). 
1.3.1 Local operations 

An operation is called local if under its action the subsystems evolve inde- 
pendently from each other. In terms of operator sums this is expressed (for 
bipartite systems, with a straightforward generalisation for the multipartite 

case) as 

8ioc(g) = E E i ® Fj g E\ ® F}, with £ E\Ei ® F}F 3 = 1 H ^ H2 . (17) 
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Local unitary evolutions Ui oc = IA\®U.2 are just special cases of general local 
operations. Since both subsystems evolve independently from each other, pos- 
sibly preexisting correlations remain unaffected. A product state will remain 
a product state, 

£ioc(Qi ® Q2) = (E EiQiEf) ® (E F i& F i) , (18) 

i i 

and any separable state will remain separable under local operations: 

Sioc (E Pi A ® ei) = E ft (E ^ei^) ® (E • ( 19 ) 

Therefore, starting from a separable state no correlations - neither classical 
nor quantum - can be created by local operations alone. 



1.3.2 Global operations 

If two subsystems are interacting with each other, their evolution will in gen- 
eral not derive from purely local operations. Any operation that is not local 
is called global. Under this type of operations all kinds of correlations can in- 
crease, as well as decrease. Therefore, entangled states can be created from 
initially separable states and vice versa. The most prominent and natural way 
of creating entangled states is a global unitary evolution due to an interaction 
between subsystems. 



1.3.3 Local operations and classical communication (LOCC) 

A prominent subclass of global operations are local operations and classical 
communication (LOCC). They comprise general local operations, and also al- 
low for classical correlations between them. The idea behind it is to allow 
arbitrary local operations and, in addition, to admit all classical means to 
correlate their application. Hence, parties having access to different subsys- 
tems can use means of classical communication to exchange information about 
their locally performed operations and the respective outcomes, and, subse- 
quently, apply some further local operations conditioned on the communicated 
information. 
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In terms of operator sums this can be expressed as 4 

£(Q) = Y,E i ®F i Q E\ <g> F}, with Yl E \Ei ® F}Fi = l Hl ® H2 . (20) 

i i 

In contrast to Eq. (17), only a single sum is involved in the description of 
LOCC operations. This is a manifestation of the correlated application of the 
respective operations on the subsystems: if the operator is applied to the 
first subsystem, the operator Fj is applied to the second subsystem. 

LOCC operations can be used to create classical correlations between subsys- 
tems. In general, a product state will not remain a direct product under the 
action of an LOCC operation: 

£locc(£i ® Q2) = E (EiQiEt) ® {F lQ2 F}) =J2 Pi Q\ ® Q\ , (21) 

i i 

with e j = E iQl Ej /Tr (E iQl Ej), q\ = F q 2 F} /Tr (f^ 2 F/), andp, = Tr (^£?i4) Tr (f^ 2 F/). 
Thus, classical probabilistic correlations can change under the action of LOCC 
operations. However, any separable state will always remain separable un- 
der LOCC operations. Accordingly, entangled states cannot be created with 
LOCC operations. 

1.4 Entanglement monotones 

Since we have argued that entanglement cannot be created using LOCC op- 
erations, our discussion at the beginning of Section 1.3 suggests to consider 
quantities that do not increase precisely under LOCC operations to quantify 
entanglement. Any scalar valued function that satisfies this criterion is called 
an entanglement monotone [25,41]. 

1.4-1 Pure bipartite states 

For pure bipartite states there exists a simple criterion that allows for the 
characterisation of entanglement monotones. It was shown that a state |S) 
can be prepared starting from a second state |$) and using only LOCC, if and 

— * — * 

only if the vector A s of Schmidt coefficients of |S) majorises A$ [47] 

A$ -< A s . (22) 

4 Strictly speaking, Eq. (20) characterises separable operations that include LOCC- 
operations. Though not every separable operation is LOCC. The exact definition 
of LOCC reads Specie) = Eijk... ■ • • (1 ® F ijk )(E i3 ® 1)(1 ® F t ) e {\ ® F^E^ ® 
l)(4 fe ®l)...[46] 
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Majorisation means that the components and [A$]j of both vectors, listed 
in nonincreasing order, satisfy X^=i[As]i > Si=i[A$]i, for 1 < j < d, with 
equality when j = d (due to normalisation, see Eq. (8)). Since the Schmidt 
vector A m , with equal components 1/d as introduced in Section 1.2.1, is ma- 
jorised by any vector A, any bipartite state can be prepared with LOCC start- 
ing out from a state \^ m ) with Schmidt vector \ m . This justifies calling |\P m ) 
'maximally entangled'. 

Since entanglement cannot increase under LOCC operations, any monotone 
M has to satisfy 

A<($) > M(E) , for A$ -< A s . (23) 

This condition is known as Schur concavity. It is satisfied if and only if [48] Ai 
given as a function of the Schmidt coefficients is invariant under permutations 
of any two arguments and satisfies 

Due to the above-mentioned invariance, there is nothing peculiar about the 
first two components of A - if Eq. (24) holds true for Ai and A 2 , it is satisfied 
for any two components of A. 

The above characterisation allows to derive several entanglement monotones 
for pure states. Very useful quantities in this context are the reduced density 
matrices, Qi or g 2 , obtained by tracing over one subsystem 

ei = Tr 2 \*)(*\, Q 2 = Tr 1 \V)(V\ . (25) 

The basic idea is that the degree of mixing of a reduced density matrix is 
directly related to the amount of entanglement of the pure state \^}. Any 
function g(g r ) of a reduced density matrix g r that is 

- invariant under unitary transformations, g(g r ) = g(Ug r W), and 

- concave, g(g r ) > Xg(g a ) + (1 — A)g(f)g), for any < A < 1, and states g a 
and gp such that g r = \g a + (1 — A)^, 

is Schur concave, and therefore provides an entanglement monotone = 
g(g r ) [25]. The most prominent choice of g is the von Neumann entropy 

S(g r ) = — Tr£ r In g r (26) 

of the reduced density matrix, often simply called the entanglement E("&) = 
S(g r ) of the pure state \^}. 

Note that due to the invariance of g under unitary transformations, g can only 
be a function of unitary invariants, hence, of the spectrum of g r . Accordingly, 
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it is not necessary to distinguish between Qi and g 2 , since they have the same 
non-vanishing eigenvalues. If both subsystems have the same dimensions, the 
spectrum of £1 equals that of g 2 - If the dimensions are not equal the reduced 
density matrix of the larger subsystem has some additional vanishing eigen- 
values. That is why one often does not distinguish between £1 and g 2 , but 
rather expresses M(^) as 

M(V) = g{Q r ) , with Qr = T^|*)(tt| , (27) 

'the' reduced density matrix, where the partial trace Tr p does not specify 
explicitly which subsystem is traced out. In subsection 2.3.1 we will discuss 
a situation where the proper choice of the subsystem over which the trace is 
performed is not completely arbitrary. 

Finally, it should be kept in mind that a single monotone is in general in- 
sufficient to completely characterise the quantum correlations contained in a 
given pure state. For such a characterisation, knowledge of all - i.e. d — 1 
(with d = min(dim(7Yi), dim(7Y2))) independent Schmidt coefficients is re- 
quired. For pure states this does not represent a serious problem but already 
indicates that, with increasing dimension, the complete characterisation of 
arbitrary mixed states will define a task of rapidly increasing complexity. Fur- 
ther down in this review, when dealing with higher dimensional mixed states, 
we will therefore have to specify which specific type of correlations we want 
to scrutinize, rather than to give an exhaustive description of all correlations 
inscribed into a given state. 

1.4-2 Mixed states 

For pure bipartite states it is rather simple to find some entanglement mono- 
tones - any unitarily invariant, concave function of the reduced density matrix 
defines one. This is due to the fact that there are no classical probabilistic 
correlations contained in pure states. For mixed states the situation is much 
more involved, because there are both classical and quantum correlations that 
have to be discriminated against each other by an entanglement monotone. 
It is by no means obvious to devise a unique generalisation of a pure state 
monotone M(^) to a mixed state monotone M.{q), such that 

- A4(g) reduces to the original pure state definition when applied to pure 
states, and 

- M.{q) is an entanglement monotone, i.e. non-increasing under LOCC. 

We will here follow one particular generalisation that applies to any pure state 
monotone [25,41], and therefore is the most commonly used one. It can be 
easily formulated, but poses severe problems when it comes to its quantitative 
evaluation. Any mixed state can be expressed as a convex sum of pure states: 
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g = J2 Pi I^X^I • (28) 

i 

On a first glance, it might appear as a self-evident generalisation to sum 
up the entanglement assigned by a certain monotone to the pure states in 
Eq. (28), weighted by the prefactors pi. Unfortunately, the decomposition into 
pure states is not unique, and different decompositions in general lead to 
different values for a given entanglement monotone. The proper, unambiguous 
generalisation of a pure state monotone, that we will also use in the following, 
therefore uses the infimum over all decompositions into pure states - the so- 
called convex roof [49] 

M(g) = inf V Pi M(*i) , with Pi >0, and g = I>i|tfi)<*i| ■ (29) 

i i 

An explicit evaluation of this quantity for a specific state - one has to find 
the infimum over all possible decompositions into pure states - implies a high 
dimensional optimisation problem - in general a very hard computational task. 

To ease this enterprise, it is convenient to make use of the following character- 
isation of all ensembles of pure states which represent a certain mixed state. 
Using sub-normalised states 

\A) = VpI (30) 

allows to reduce the number of involved quantities. Since the Pi are positive, 
one has — Pi\^i){^i\. Assume one ensemble {iV'i)} is known such that 

Q — Hi iV'iKV'il " e -9- the eigensystem of g . New ensembles defined as 

10.) = E with £ VlV l3 = 5 jk , (31) 

3 i 

represent the same mixed state g = J2i I^X^I = Hi and any ensem- 

ble representing g can be constructed in this way [50,51]. In the following, any 
matrix V satisfying Eq. (31) will be referred to as left unitary. The number, 
cardinality, of ensemble members of the decomposition of g (i.e., the length 
of the index set of % in Eq. (31)) is not fixed by the rank, i.e., the number of 
nonvanishing eigenvalues of the density matrix which represents the consid- 
ered state. There is no a priori maximum cardinality, though it is sufficient 
to consider ensembles with cardinality not larger than the square of the con- 
sidered state's rank [52]. However, there is no evidence that it is necessary to 
employ ensembles of this maximum cardinality, and, in particular, there is no 
proof that the infimum in Eq. (29) cannot be found with smaller ensembles. 
Nonetheless, without a sharper bound on the length of the decomposition, we 

2 

need to find the optimal left-unitary matrix V G C n xra , which implies an op- 
timisation procedure of dimension ~ n 3 to compute the entanglement A4(g) 
of a given state g of rank n. Since there is no simple parametrisation of ar- 
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bitrary left-unitary matrices, the constraint J2i VkiYij = $jk even complicates 
numerical implementation. 



1.5 Entanglement measures 



Entanglement monotones that satisfy some additional axioms are called en- 
tanglement measures E. So far, however, there is no uniquely accepted list of 
axioms, hence there is no commonly accepted distinction between monotones 
and measures [53,54]. We do not attribute too much relevance to this question 
of terminology, and just present here a list of potential axioms: 



- E(g) vanishes exactly for separable states. 

- additivity: the entanglement of several copies of a state adds up to n times 
the entanglement of a single copy, E(g® n ) = nE(g). 

- subadditivity: the entanglement of two states is not larger than the sum of 
the entanglement of both individual states, E(g <g> g') < E(g) + E(g'). 

- convexity: E(Xg + (1 - X)g') < \E(g) + (1 - X)E(g'), for < A < 1. 



Some authors additionally require that an entanglement measure has to be 
invariant under local unitary transformations. However, this is already im- 
plied by monotonicity under LOCC (Section 1.4), since monotonicity implies 
invariance under transformations that are invertible within the class of LOCC 
operations [55]. Local unitaries and their inverses are LOCC operations. Thus, 
any entanglement monotone and measure has to be non-increasing under both 
the former and the latter. Since non-increasing behaviour under the latter im- 
plies non-decreasing behaviour under the former and vice versa, any monotone 
or measure has to be invariant under local unitaries. 



There are attempts to find a distinct set of axioms that leads to a unique 
measure [56]. On the other hand, it is sometimes necessary to relax some of 
the above listed constraints, in order to find a measure that is computable. 
For example, negativity [26] has become a commonly used quantity although it 
vanishes for a class of entangled states [26] . Though, compared to other mea- 
sures it has the major advantage that it can be computed straightforwardly. 
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2 Concurrence 



2.1 Two-level systems 



Concurrence was originally introduced as an auxiliary quantity used to calcu- 
late the entanglement of formation of 2 x 2 systems. However, concurrence can 
also be considered as an independent entanglement measure [27]. The origi- 
nal definition of concurrence [41,57] for bipartite two- level systems is given in 
terms of a special basis 



| ei > = |$+>, |e 2 > =*|$-> , |e 3 )=^+), |e 4 > = |*~>, 



(32) 



where |$ ± ) = (|00) ± |11>) /V2 and = (|01) ± |10))/V2 are the Bell 

states [41]. Using this particular basis, the concurrence c of a pure state 
is defined as 



cm 



Writing this definition more explicitly, c(\I/) 
after summation, with the alternative formulation [27] 



(33) 

Ei(**|eJ')(e i |*) , one ends up, 



cm 



{V*\<Ty ® <Ty\V) 



(34) 



where o~ y is the second Pauli matrix, and \^/*) is the complex conjugate of |\&) 
with the conjugation performed in the standard (real) basis { 1 00), |01), |10), |11)}. 
Since a scalar product implies a complex conjugation anyway, the second con- 
jugation cancels the first one such that (\&*| = is the transpose and 
not the adjoint of = Eq. (34) is the most commonly used for- 
mulation and is often considered as the definition of concurrence rather than 
Eq. (33). 

The concurrence of mixed states is given by the corresponding convex roof, 
alike Eq. (29): 

c(g) = inf Pi c (^i) , with Pi > , and g = ^2pi\%){%\ ■ (35) 

Concurrence of a pure subnormalised state \ip) (see Eq. (30)) can be expressed 
as c(tp) = \f c (ip,ip)\, in terms of the function 



31 i 



(36) 



that is linear in both arguments. These linearity properties, together with the 
parametrisation of all decompositions of g into pure states given in Eq. (31), 
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allows to write Eq. (35) as 



V T ~ 






ki 



(37) 



The quantities f c (ifjj,ipk) can be understood as elements of a complex 
symmetric matrix r. Hence one can use the compact matrix notation 



c(£) = inf]T [VrV 1 

i 

The infimum of this quantity is known [27] to be given by 

c(q) = max j«Si - °| , 



(38) 



(39) 



where the Si are the singular values of r, in decreasing order. They can be ob- 
tained as the square roots of the eigenvalues of the positive hermitian matrix 



tt"I" . 



Since we will refer to infima of expressions similar to that given in Eq. (38) 
several times later on, we discuss the derivation of this infimum in some detail. 
Though, we do not follow here the original derivation [27], but rather present 
a generalisation [49] valid in arbitrary dimensions - which we will need for 
later reference when considering subsystems with more than two levels. Thus, 
the following considerations do not only apply to the hitherto discussed 2x2 
case, but also to systems of arbitrary finite dimension. 



Any complex matrix M G C niX ™ 2 can be diagonalised [35] as 

M = UiDU r , (40) 

where U\ G (D" lXd and U r G (C dxn ' 2 , with d = min(ni, 722), are, respectively, left- 
and right unitary, i.e. U}Ui = 1 and U r U}. = 1, and D G <C dxd is a diagonal 
matrix with real and positive diagonal elements, referred to as singular values 
of M. Moreover, C/j and U r can always be chosen such that the singular values 
are arranged in decreasing order along the diagonal. 

Applying this to a square, complex symmetric matrix r, one concludes that r 
can be diagonalised with a unitary transformation U as 

f/rf/ T = diag[ l S 1 ,...,5 n ] =:r dia S. (41) 

Given this diagonal representation, one defines a transformation Vh with the 
help of n x 2 k Hadamard Matrices H [58] with 2 k > n. A square 2 k x 2 k 
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separable state entangled state 




Si 3? Si 3? 

Fig. 1. Complexified singular values Si and —Sje 1 ^ 3 (j > 1) - see Eq. (42) - plotted 
in the complex plane. For a separable state (left) one can always find appropriate 
phase factors such that all terms add up to 0. However, for an entangled state (right) 
this is not possible. The minimum of Eq. (42) can be obtained for e i(p j = 1 for all j. 



Hadamard matrix can be constructed for each k and, by definition, has its 
columns given by mutually orthogonal real vectors, with the same absolute 
value l/(2 k l 2 ) of all elements. We shall denote by H (and call it also a 
Hadamard matrix) a rectangular n x 2 k matrix obtained from the original 
square one by keeping only n rows. Due to the rows' orthogonality, H is left 
unitary, WH = t. 



The transformation matrix Vh is a modification of a Hadamard matrix. Namely 
each column j but the first one is multiplied with a phase factor ie ttpj / 2 . The 
latter does not affect left- unit arity, VjjVh = 1. However, since enters 
Eq. (38) instead of V#, the phase factors are indeed important. Carrying out 
the transformation, one obtains 



t/ ^diag T /T 



1 

2^ 



(Si-^e^) , 



(42) 



where one needs not care about the non-diagonal entries, since only diagonal 
elements are summed up in the end. 

So far we found a (left- unitary) transformation VhU such that 



E 



[v h UtU t V£ 



(43) 



Now one has to distinguish two cases. In the first case one has Si < X^>i Sj- 
In that case one can always find phases ipj such that Si — J2j>i Sje lLP: > = 0, as 
depicted in Fig. 1. In the second case, where «Si > J2j>i Sj, the optimal choice 
for all phases (such that Si — J2j>i Sje llp i be minimal) is (fj = 0, and one gets 
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Si — J2j>i$j- Altogether, we found a transformation such that 



£ [v H UrU T V£ 



max < Si — ^2 



(44) 



It is now easy to show that there is no left unitary transformation lead- 
ing to a smaller result. We can restrict ourselves to the second case above, 
<~>i > Sj>i since is the smallest possible value of a non-negative quantity 
anyway. To do so, we start with the diagonal form r^^ a S of r, and transform 
it with Up = diag(l, i, . . . , i) such that we obtain 



f dia S = U p r di ^Uj = dia^ 



Si , — 02 , • • • , — S n 



Now one has 



EE* 



~diag 



ij 'jk 



i jk 



> 



E(K?5i-E^ 3 

i j>i 



(45) 



(46) 



One can always choose V in such a way that the Vn are real for any %. Then 
one has J2iV£ = 1, and |Z)jV^| < 1. Therefore: 



EE* 



,diag 



(47) 



i>l 



Thus we have shown that Eq. (39) really is the infimum corresponding to 
Eq. (38). 

From this exact expression for the concurrence for mixed states, one can also 
deduce the entanglement of formation. As we have seen in Section 1.4, the 
entanglement E(ty) of a pure state |\&) can be quantified by the von Neu- 
mann entropy S(g) of the reduced density matrix, Eq. (26). Entanglement of 
formation of a mixed state then follows as the convex roof [41] 



E(q) 



{p 



inf V Pi E(Vi) , with Pi > , s.t. q = ]T Pt 



(48) 



In arbitrary dimensions the underlying optimisation problem is unsolved - 
apart from a few known solutions for particular states [59]. Only for 2 x 
2 systems an algebraic solution is known for general states. For these low 
dimensional systems, the entanglement of a pure state can be expressed as a 
function of its concurrence, E(ty) = £(c(\I/)), with [27] 

£(c) = - £ i(l + «v / T^))log 2 i(l + av / r^)) . (49) 

a=— 1,1 



The function 5(c) is monotonically increasing and convex. Thus, the entan- 
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glement of formation can be estimated (with (29)) as 

E{q) = inf ]T Pl 8(a) > inf sfe Pl c t ) = £(c(g)) . (50) 

i i 

Consequently, concurrence provides a lower bound for the entanglement of for- 
mation. In general, the decomposition that provides the infimum is not unique. 
In the present case of 2 x 2 systems, the manifold on which the infimum is 
adopted always contains a set of pure states all of which have the same concur- 
rence [27]. For these special decompositions one has J2iPi£(°i) — ^(j2iPi c ij ■ 
Thus, equality holds in the above inequality. Therefore, in 2 x 2 systems, one 
can always express entanglement of formation in terms of concurrence. Since 
an algebraic expression, Eq. (39), for concurrence is available for arbitrary 2x2 
states, also entanglement of formation can be computed purely algebraically. 

2.2 Higher dimensional systems 

In the following Section we will focus on the concurrence of systems of arbitrary 
dimension. Entanglement of formation is more appealing than concurrence, 
because it is believed [60,61,62] - though not proven for general states - to 
be additive, whereas it is evident that concurrence is not additive. However, 
concurrence satisfies several algebraic properties that provide a basis for good 
approximations, while it is unknown whether entanglement of formation can 
be evaluated as efficiently Thus, whereas entanglement of formation may be 
more appealing from a rather fundamental point of view, concurrence is more 
appealing for pragmatic reasons - it allows for the efficient description of states 
even in high dimensional systems, as we will see in the following. 

The quantitative estimation of the concurrence of a mixed state is often 
achieved by numerical means [63,64,65] which essentially solve a high dimen- 
sional optimisation problem when searching for the minimum that defines the 
convex roof, Eq. (29). However, such an approach can only provide an upper 
bound of concurrence rather than its actual value, since a numerical optimi- 
sation procedure can never guarantee convergence into the global rather than 
into a local minimum. Hence, besides the efficient numerical implementation 
of optimisation procedures, it is most desirable to derive a lower bound of 
concurrence that can be possibly evaluated in a purely algebraic way, since 
the numerical effort increases rapidly with the dimension of the underlying 
Hilbert space. In the present Section, we will derive an approach to charac- 
terise concurrence of general mixed states in arbitrary, finite dimensions. In 
particular, we will 

- not only provide a framework for an efficient numerical implementation to 
compute an upper bound of concurrence in Section 2.3.4, but also 
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- propose a generalisation of concurrence for multipartite systems in Sec- 
tion 2.3.5. 

- Moreover, we will formulate lower bounds of concurrence in Section 3, some 
of which can be computed purely algebraically. 

- Finally, in Section 3.2 we will derive an approximation of concurrence that 
is valid for most states describing current experiments, and which can also 
be evaluated purely algebraically. 

Note that, when describing higher dimensional systems, our choice of concur- 
rence as a single scalar quantity can never explore all correlations inscribed in 
an arbitrary given state (see also the discussion at the end of Section 1.4.1). 
Therefore, the definition we shall elaborate here will be constructed such as 
to allow to target at different, specific types of correlations, in a possibly mul- 
tipartite, higher dimensional quantum system. 

Both, upper and lower bounds of these concurrences, will allow to confine their 
actual values to finite intervals, providing reliable information about arbitrary 
states. Specializing our approach, in Section 3.2, to typical experimental re- 
quirements, will relax the demand for a completely general treatment, and will 
be rewarded by a dramatic speed-up of actual numerical evaluations, through 
a very efficient and easily implemented estimate of concurrence. 

Since we have neither an a priori estimate of the tightness of our bounds, 
nor one for the range of validity of our approximation, we will compare our 
estimations from below with the corresponding upper bound in Section 3.3. 
For this purpose, we will use random states, as well as states under scrutiny 
in real experiments [9]. 

To start with, we have to realize that the definition of concurrence given in 
Section 2 only applies to two-level systems. Since Bell states used in the orig- 
inal definition in Eq. (33), and the spin flip operation used in Eq. (34) do not 
have unique generalisations to higher dimensions, there is no straightforward 
generalisation of concurrence to higher dimensions. So far, two inequivalent 
generalisations for systems comprising more than two levels have been for- 
mulated, which both coincide with the original one, if restricted to two-level 
systems. 

2.2.1 Q-concurrence 

One possible generalisation is Q-concurrence [49]. The complex conjugation 
that appears in the bra in Eq. (34) can be perceived as an anti-linear operation. 
G-concurrence is based on anti-linear operators, where anti-linearity of an 
operator 6 is defined by the property 

e(ai|*i) + a 2 |* 2 )) = aje|*i) + a;0|^ 2 > ■ (51) 
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An anti-linear operator that is unitary _1 = ©^ and an involution = _1 
is called a conjugation. In terms of such a conjugation 0, one can define 0- 
concurrence ce(V0 °f a pure ( n °t necessarily normalised) state \ip) as 

c©(V) = |<V|ej^>| . (52) 

Of course, c@(ip) does not only depend on but also on the choice of 0. Thus 
©-concurrence is not a single, uniquely defined quantity, but rather a family 
of quantities, depending on the choice of 0. In systems larger than two- level 
systems, no conjugation is known such that ce(V') vanishes for all separable 
states and is strictly larger than zero for all entangled states. However, in the 
case of two-level systems there is one. For = o~ y <g> a y [49,66], with the 
second Pauli matrix a y and C* the complex conjugation in the standard basis 
(defined after Eq. (34)), c@(ip) is non-vanishing if and only if \ip) is entangled. 
For this special choice, ©-concurrence coincides with regular concurrence as 
defined in Eq. (34). 

©-concurrence can easily be extended to mixed states using the concept of 
convex roofs. Given a complex symmetric matrix t@ with elements 

[re]jk = (iij\Q\ipk) , (53) 

one easily finds for the ©-concurrence cq(q) = inf J^iPiCe^i) of a mixed state 
q: 

(54) 



c e (£) = inf£ [Vr e V 



As discussed in Section 2, Eq. (39), the infimum can be expressed as 

c (f?) = max| l S 1 e -£«Sf, oj , (55) 

with the singular values Sf of r in decreasing order. Thus, ©-concurrence 
can be easily evaluated for arbitrary mixed states. However, since, apart from 
= o- y ®o~ y C* in 2x2 systems, no conjugation is known that is positive exactly 
for entangled states, it has the disadvantageous property that c vanishes for 
some entangled states. 



2.2.2 I- concurrence 

Alternatively, I-concurrence [67] is defined in terms of operators I\ and I2 
acting on 13 (Hi) and B(H 2 ) as 



c/(#) = /^M^)Wf> • (56) 
The operators U are required to satisfy the following properties [67] 
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a) IiH = (htiy (i = 1,2), for all hermitian operators H, which ensures that 
/-concurrence is real. 

b) [Ii,U] — (i — 1,2) for all unitary U, which ensures that /-concurrence is 
invariant under local unitary transformations. 

c) (^|(/i <g> / 2 |*)(*|)|*) > 0, for all states |^), where equality holds if and 
only if |\&) is separable . 

Up to scaling, there is a unique operator satisfying these requirements [67], 
namely 

= % - \<Pi){<Pi\ ,\<Pi)eHi, i = l,2, (57) 
which maps \ifii)(ipi\ onto its orthogonal space. Thus - in contrast to 0- 
concurrence - /-concurrence is a quantity that is uniquely defined up to a mul- 
tiplicative constant. With Eq. (57), /-concurrence cifi?) of a pure state \^f) can 
also be expressed in terms of reduced density matrices. With Ii <g> I 2 \ty){$ \ = 
1 — Qi <g) 1 2 — li <8) Q2 + |^ r )(^ r |, one easily obtains 

Cl (g) = v/2-Tr^-Tr^ . (58) 

As argued right before Eq. (27), the last two terms are equal, such that there 
is no need to explicitly distinguish between the two reduced density matrices. 
It therefore became a widespread convention to define concurrence using only 
one of the two reduced density matrices 

c(tt) = ^2(1-1^2) , (59) 

where g r can be either one. 

If we now use the Schmidt form, Eq. (8), of an arbitrary pure state |\&), its 
/-concurrence reads 

c,(*)= /2]T SiSj, (60) 

V <#j 

and it is easily verified that /-concurrence coincides with the original definition 
given in Eq. (34), for two- level systems. Note that /-concurrence cannot exceed 
a given maximum value, 

cj^) < v/2(l - 1/d) , (61) 
with d the dimension of the smallest subsystem. 

2.3 Representation in product spaces 

To start with, we would like to represent the definition of concurrence, Eq. (56), 
in a different form. Whereas Eq. (56) allows for a nice interpretation of the 
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first second 
subsystem subsystem 



first copy 




second copy 



o 




p(l) 



p(2) 



Fig. 2. An expectation value (^|a|^) with respect to a single copy of a given state 
|*) can not discriminate separable from entangled states. However, this is possible, 
if expectation values with respect to two copies are used. The grey boxes symbolise 
the two copies of the state |^), with first subsystem (black) and second subsystem 
(grey). This scheme defines the structure of an operator acting on two copies of the 
single subsystems, such that the corresponding expectation values vanish exactly if 
the state |^) is separable. 

entanglement of pure states |^) in terms of the degree of mixing of the re- 
duced density matrices, it is not very suitable for the evaluation of the convex 
roof (29), if we are dealing with mixed states. What we are looking for instead, 
is a linear operator a such that all expectation values (^lal^) with respect to 
separable pure states vanish, whereas they are strictly positive for entangled 
pure states. However, such an operator does not exist. In the case of two- level 
systems, one could circumvent this problem by using ($f*\ instead of (see 
Eq. (34)). However, this trick does not work properly in higher dimensional 
systems, as discussed in 2.2.1. 

In order to obtain a quantity that unambiguously detects all entangled states 
in arbitrary dimensions, we can follow a different way. We consider expectation 
values with respect to two copies of the pure state under investigation, 



where A is acting on 7i <g> H, i.e on Hi ® H2 ® Hi <S> 7^2- Of course, one needs 
to require that Eq. (62) vanishes for any separable state = \ip) <g> \(f>). 
The simplest possibility could be an operator A that decomposes into a part 
acting only on the copies of the first subsystem and a second part acting on 
the copies of the second subsystem. Not necessarily the unique, but a good 
choice for A are projectors P_ ^ (k = 1,2) onto antisymmetric subspaces Ji k _ 
of the space Tik <8> Ti-k- They contain all states that acquire a phase shift of n 
under the exchange of the two copies of Tik- Thus, any state | ■?/>_) G 7i k _ can 
be expressed as \ipJ) = Vfy'Ou) — where the states {\i)} and {\j)} 
form an arbitrary basis of Hk- 




(62) 
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Indeed, since the two-fold copy \ip) ® \ip) of a state is a symmetric object - 
it remains invariant under an exchange of the copies - the expectation value 
(ip\ ® (ip\P^\ip} ® \ip) vanishes for any state \ip) G Hi, and the same holds 
true for the analogous expression for H 2 . 

Now, one can define A := 4P a in Eq. (62), with the projector P a onto the 
space spanned by the states in Hi®H 2 ®Hi®H 2 that are antisymmetric both 
with respect to an exchange of the two copies of Hi, and with respect to an 
exchange of the two copies of H 2 . The concurrence in Eq. (62) then necessarily 
vanishes for separable states. The prefactor 4 is just a normalisation chosen 
such that the concurrence ranges from to 1 for two-level systems. With this 
normalisation, the present definition is indeed equivalent to Eq. (59). 

One may find it easier to interprete A in terms of the two copies of the single 
subspaces Hi and H 2 . For this purpose, let's identify Hi®H 2 ®Hi® H 2 and 
Hi®Hi®H 2 ®H 2 . Then - as illustrated in Fig. 2 - A is just the tensor product 
of the two projectors onto the two involved anti-symmetric subspaces: 

A = AP { } ] <g> Pl 2) . (63) 

If a state \%) is separable, the expectation value (^\ <8> ® fac- 

torises into the product of the analogous expressions corresponding to the 
single subsystems 

(tf | <g> (V\A\V) <g) I*) = 4{<p\ <g> {(p\P^\(p) ® \(p){</>\ ® ® 10) , (64) 

and vanishes because both, \ip) <8> \<p) and \<p) <g> \<p), are symmetric. However, 
for an entangled state |\&) = J2i V~K\fii)i <S> \£i)2 - for convenience represented 
in its Schmidt decomposition, Eq. (8) - the two-fold copy <g) |^) is not 
symmetric under exchange of the copies of any of the subsystems. Thus, <g> 
|\&) necessarily has some anti-symmetric part, and the expectation value of A 
with respect to this two- fold copy of \^/) is strictly positive. 

Now, with this definition of A, one can also rephrase the concurrence of a 
mixed state g in terms of subnormalised states, Eq. (30), as 

c (^) = r f n f E v^il l^i) ® ■ (65) 

If one makes use of the prescription (31) to characterise all decompositions 
of £> into pure states, it reveals useful to define a tensor A that contains the 
elements of A evaluated with the subnormalised states 

A 1 ™ = ® (0 m | A <g> |0 fc ) . (66) 
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One can now easily rewrite c(g) in the closed expression 



c(q) = mfJ2]/[V® V A , (67) 

i 

where the optimisation is to be performed over all left unitaries V. 
2.3.1 Symmetries of A 

We required A to be anti-symmetric with respect to the exchange of the two 
copies of Hi as well as of H 2 - Though, this choice is not unique. It would have 
been sufficient to require only one of these symmetries. Here, we will briefly 
sketch the consequences of specific choice of A. In Section 2.3.2 and later on 
in Section 3, it will become apparent why it is of importance to require both 
symmetries. 

With our above definition, Eq. (63), of A, the elements of A 
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TV, 



\<f>k)((f>r. 



Tr 9 



Th Tn \<t> k ){<t> m \ +Tr |^)(<A 



Tr 



l0JfeX0r. 



(68) 



contain partial traces over both subsystems, in a balanced way, what would 
not have been the case for other choices. E.g., for a projector A onto the space 
spanned by all the states that are antisymmetric with respect only to the 
exchange of TC\ , the analogous expression reads 



Af k =2 



Tr 



\4>j){<t>l\<t>k){<t>r, 



Tr x 



Tr 5 



\<t>j){<t>l\ Tr 2 |</>fcX</> 



(69) 



Whereas the symmetric treatment of both subsystems does not have any prac- 
tical consequences, and is of rather aesthetical character, there is a much more 
crucial property: A is invariant under exchanges of the co- or contra-variant 

A" 1 }!. This symmetry will turn out to be the crucial 



indices, A 1 -]! 



A lm 



ingredient required for the approximations to be discussed in Section 3. 



2.3.2 Two-level systems 

In the case of two-level systems, there is only one ant i- symmetric state, namely 
1 01 } — 1 10). Therefore, the projectors onto the anti-symmetric subspaces H}_ 
and 1~L 2 _ have only one non-vanishing eigenvalue. This, in turn, also holds true 
for A which reads A = |x)(x|, with 

|x) = |0011) - |0110) - |1001) + |1100) . (70) 
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Consequently, A can be expanded in terms of a single matrix r as 

•^t = rLrjk , with Tjk = ( X \<t>j) ® \<Pk) ■ (71) 

Due to the symmetry of A discussed above in Section 2.3.1, r is indeed sym- 
metric, i.e. satisfies a crucial precondition for the generalisation to mixed 
states. Expressed in terms of r, Eq. (67) now simplifies to 



<Q) = i£ f £ y/[V^Vn u [VtVT\ u = inf I [VtV t ] J . (72) 

i i 

Due to the symmetry of r, the infimum is exactly given in terms of the singular 
values of r, as discussed in the context of Eq. (38). 

The relation between the present and the original approach to concurrence [27,57] 
becomes apparent with the observation that the matrix r defined here coin- 
cides with that introduced in Eq. (38) 

T jk = (X\<l>j) ® \<f>k) = (<t>jWy ® <Ty\<l>k) , (73) 

in the 2x2 case. Thus, the original approach to concurrence is naturally 
embedded in the present, more general framework. 



2.3.3 Higher dimensional systems 

The formalism in terms of projectors onto anti-symmetric subspaces can also 
be used to formulate a generalisation of concurrence to systems of arbitrary di- 
mensions. For an n\ x ^-dimensional system, Ah&sm = ni(ni — 1)^2(^2 — 1)/4 
non- vanishing eigenvalues - ni(ni — l)/2 corresponding to the antisymmetric 
subspace of Hi <8>Hi, and analogously for H 2 <E> H 2 - Thus, A cannot be ex- 
pressed with the help of a single eigenvector anymore, but with a finite sum 

m 

A=Y,\Xa)(x«\, (74) 
0=1 

which contains ni{n\ — l)n 2 (n 2 — l)/4 non- vanishing terms. Due to the sym- 
metries A l j% = A l k ™ = Ajk which we already observed in Section 2.3.1, all T a 
with the elements 

Tf k = (Xa\<f>j) ® \<f>k) ,a = l,...,m = m(ni-l)n 2 (n2-l)/4, (75) 

are symmetric. Furthermore, the eigenvectors of A still carry an undetermined 
phase factor e tlfa . Whereas these free phases usually do not matter, they pro- 
vide an additional freedom which we shall exploit in Section 3. Therefore, we 
explicitly account for the free phases <p a , and Eq. (67) consequently reads 

<Q)=tyEJE\\y^* v ~ vT \*\ 2 - ( 76 ) 

i y a 
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There are two crucial differences as compared to Eq. (72) that hamper finding 
the exact infimum: First, the square and the square root in Eq. (76) lead to 
non-linear expressions in the T a , and, second, the already mentioned fact that 
the several distinct symmetric matrices T a cannot, in general, be diagonalised 
simultaneously. 

One of the earliest generalisations of concurrence to higher dimensional sys- 
tems that does not lead to the non-linear behaviour, as it appears in Eq. (76), 
is the concurrence vector. Although the original definition [64] is slightly dif- 
ferent, we will describe it here in terms of the projectors of Eq. (63), in order 
to highlight the similarities with our own approach. Each \\ a ) inherits a neg- 
ative pre-factor under the exchange of the two copies of either Tii or Ti,2- The 
product |^ s ) <S> \^ s ) °f an arbitrary separable state is invariant under such 
transformations. Since the overlap 

C Q (tt) = ( Xa m ® |*> (77) 

is invariant under the exchange of the two copies, it necessarily needs to vanish 
for any \x a ) and any separable state = \^/ s ) <8> \^ s )- In contrast, whenever 
C a (^>) is positive, the state \%>) necessarily needs to be entangled. The inverse 
implication is a bit more involved - |\&) is separable, if C a (^>) = 0, for a = 
1, . . . ,ni(rai - l)n 2 (n 2 - l)/4. 

The generalisation of the above for mixed states is now straightforward and 
analogous to the case of two-level systems discussed in Section 2.3.2. A given 
state g is separable if and only if there is a left-unitary transformation V such 
that all elements 

(78) 

of the concurrence vector C vanish. However, this - necessary and sufficient 

- separability criterion is, in general, difficult to evaluate - since, again, in 
general the matrices T a cannot be diagonalised simultaneously. Nevertheless, 
it establishes the basis for some operational, though only necessary separability 
criteria. It implies that a given state is entangled if the singular values 

of one matrix T a satisfy S[°^ — Z^-Sj^ > 0. Another - in general stronger 

- criterion is obtained with the help of linear combinations J2 a z aT a of all 
matrices T a , with complex pre-factors z a . For a suitably chosen set {z a }, 
the expression Si — Y,j <Sj can be significantly larger than the corresponding 
expression for a single matrix T a [64]. 

Concurrence vector is not a measure and as such does not provide an adequate 
tool for our aims, i.e. the quantitative characterisation of temporally evolving 
entanglement. Though, in section 3 we will show that - despite its non-linearity 

- Eq. (76) can be used to derive some means for such a quantification. 



c « = E [vr a v T 
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2.3.4 Gradient 



Before we focus on this, however, let's discuss how to assess concurrence nu- 
merically. As mentioned in Section 1.4.2, the cardinality of the ensemble that 
realises the infimum in Eq. (29) can exceed the rank r(g), though is bounded 
by r 2 (g) [52]. This is the cause of the appearance of a rectangular, left-unitary 
matrix V instead of a quadratic, unitary matrix U in Eq. (67). In the present 
section, however, matrices of the latter type are more convenient. Therefore, 
we will fix the cardinality of the considered ensembles. If it turns out that the 
assumed cardinality is not large enough, one can always increase it by adding 
some null-vectors to the ensemble. 

According to Eq. (67), the concurrence of a mixed state g is given by 

c(g) = inf C(U) , with C(U) = ]T yf[U®UA W ® W]$ . (79) 

i 

If one considers an infinitesimal transformation dU = 1 + ieK , with e in- 
finitesimally small, and uses the symmetry of A with respect to an exchange 
of co- and contravariant indices, this can be written as 

i — if 

C{dU) ~ V J A* + -= [K®1A -AK® l\l , (80) 

• , / A 11 

where the expansion y/a + e ~ y/a + e / of the square root function was 

employed. This allows to rephrase C(dU) as 



C(dU) = J2\fA?i + £ T, K ij G a > with G ji = i[-fi=-- f iL .(81) 




I IJ 

Since the square root is not analytic in the origin, one has to take care that 
A]] be non- vanishing for all %. This is the case if and only if there are only non- 
separable pure states in the decomposition of g, since, by Eq. (62) and 
Eq. (66), the elements A % * are the squares of the concurrences of the states 
l^i). Though, even if the denominators in Eq. (81) vanish, the numerators 
behave accordingly [68], such that Eq. (81) is indeed always well defined. 

The hermitian matrix G in Eq. (81) can be considered as a gradient, since the 
increment of C reads 

C(dU)-C(l) =e Tr [KG] , (82) 

what is just a Hilbert-Schmidt scalar product [35]. Accordingly, the direction 
of steepest descent of C is given by — G. Minima of C can therefore be found 
by repeated application of the transformation dU = exp(— ieG). However, also 
more refined methods can be used, such as the conjugate- gradient minimisa- 
tion [64,69]. In this method, the iteration is not performed exactly along the 
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current gradient G n , but it takes into account corrections that ensure that the 
previous iteration along the gradient G n is not reversed. More explicitely the 
iteration is performed along 

G n = Gi-^Gn-! , with 7 = , 

and rG n _i = e^^G^e'^- 1 . (83) 

However, there is in general no a priori information available on whether the 
solution reached by that procedure is a local or a global minimum. While one 
may start the iteration with different initial conditions parametrised by U, 
in order to get a better intuition on the effective "landscape" defining the 
optimisation problem at hand, this uncertainty persists, and the more so the 
higher the dimension of the parameter space over which the optimisation is 
carried out. 

2.3.5 Multi-partite systems 

Since lately several experimental groups [7,8,9,10,11,12,13,14] systematically 
investigate quantum correlations in multipartite systems, i.e., systems with 
more than two subsystems, a quantitative description of multipartite entan- 
glement is highly desirable. Multipartite systems keep room for distinct classes 
of quantum states: Even for a pure state |\&) of a tri-partite two- level system 
more than a single scalar quantity is required for complete characterisation of 
all inscribed quantum correlations. 

In bipartite systems, any state can be prepared using LOCC only, starting 
with a distinguished, maximally entangled state (see Section 1.4.1). This is no 
longer true in multipartite systems - where inequivalent kinds of multipartite 
entanglement exist [70]. Consider for example a Greenberger-Horne-Zeilinger 
state (GHZ-state) [22] 

|*ghz) = ^(|000) + |111)) , (84) 

and a W-state [70] 

|^ W ) = -L(|001> + |010> + |100>) . (85) 

Both states contain fundamentally different correlations, such that none of 
the two can be created from the other one by LOCC alone [70]. Thus, one 
cannot expect that a single scalar quantity can completely describe p-particle 
correlations (with p > 2). 
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Fig. 3. The concurrence of bipartite states is denned in terms of two operators, each 
of which is associated with two copies of one subsystems, as displayed in Fig. 2. This 
concept can straightforwardly be generalised to multipartite systems. For each two 
copies of one subsystems, there is one operator . 

Hence, how can we generalise the concept of concurrence to multipartite sys- 
tems? The definition of concurrence in Eq. (59) is not very suggestive for 
generalisations: While the degree of mixing of the reduced density matrix has 
an unambiguous interpretation for bipartite systems, its meaning is unclear 
in the multipartite case. Though, Eq. (62) has a rather obvious multipartite 
formulation, what implies a generalisation of concurrence for multipartite sys- 
tems that can describe at least some of the correlations we are seeking for. A 
generalisation of concurrence [71] for tri-partite two-level systems is already 
available - it characterises all tri-partite correlations of pure states. In the 
following, we present an even more general framework applicable to systems 
with an arbitrary number of subsystems, for pure and mixed states. 

Similarly to the case of bipartite concurrence, also its multipartite generalisa- 
tions can be defined in terms of an operator A acting on the tensor product of 
7i with itself. The only difference being that 7i is the tensor product of more 
than two factor spaces Tii. Any tensor product of projectors onto symmetric 
and anti-symmetric sub-spaces is positive semi-definite, and invariant under 
local unitary transformations, see Fig. 3. However, not all such operators fi- 
nally lead to tensors A with the desired invariance under exchange of co- or 
contravariant indices. This symmetry is only valid for products of an even 
number of projectors onto anti-symmetric subspaces, possibly multiplied 
with some projectors P+* onto symmetric subspaces. On the other hand, all 
expectation values of the type appearing on the right hand side of Eq. (62), 
with A a product of an odd number of projectors onto anti-symmetric sub- 
spaces, vanish anyway, for arbitrary states. 

Thus, we define iV-partite concurrence as in Eq. (62), i.e. 

c(*) = y/(V\ <8> A |¥) ® |¥) , (86) 
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in terms of a sum of direct products of projectors onto symmetric and anti- 
symmetric subspaces 



N 

A = E PM®^? , P{ Si }>°- (87) 

V{si = ±} j=l 

Here, V{sj = ±} represents all possible variations of an iV-string of the sym- 
bols + and — , and the summation is restricted to contributions with an even 
number of projectors onto anti-symmetric subspaces. 

Note that, because of the continuous parametrisation in terms of the P{ Si }, 
Eq. (87) actually defines a continuous family of multipartite concurrences. 
While the intuitive interpretation of the concurrence defined by an arbitrary 
choice of the P{ Si } so far remains an open problem, there are some specific 
choices of the P{ Si } that have immediate applications in the characterisation 
of the multipartite entangled states dealt with in Section 4. Let's take all 
prefactors p{ Si y in Eq. (87) equal, with the only exception of setting p+... + = 0. 
As in Eq. (63), there is some freedom in the normalisation. Once again we set 
P{ Sl } = 4 - just to be consistent with Eq. (63). The concurrence cjv, defined 
with these prefactors for systems with an arbitrary number iV of subsystems, 
can be expressed in terms of the reduced density matrices Qi [32] as 



c N m = 2 1 ~ ] j(2» _ 2)<*|*> 2 - £Tr ft 2 • (88) 

The index % runs over all (2 N — 2) nontrivial subsets of an iV-particle system. 
Obviously, cjy vanishes only for completely separable iV-particle states, since 
all reduced density matrices are simultaneously pure only for these states. 
What is less obvious but at least equally important is that cn allows to com- 
pare the entanglement of N- and (N — l)-partite states: for any state that 
factorises into a product of a one-particle state and a (N — l)-partite remain- 
der \<P), cjv(^) reduces to c N -i(<P). Thus, c N is perfectly suited to investigate 
the scaling properties of multipartite entanglement with the system size N, 
what we shall explore in Section 4.2.3. 

Further entanglement-properties can be addressed by other choices of the p{ Si y 
in Eq. (87): for example, bi-separability of tri-partite states can be charac- 
terised by cf 3 \ with = AP^ ® <g> P^\ and, analogously, by c^ 3 \ 
with 4 13) = 4P_ ) ® P+ ) ® P m , or by c^ 12) , with A { ^ 2) = 4Pl 1} <g> p! 2) <g> p| 3) . 
The subscript represents the number of subsystems as in Eq. (88) while the 
superscript stands for the subsets in which quantum correlations are mea- 
sured 5 . Consider |^) = |<2> 12 ) ® \ip3) as an exemplary state. The choice 

5 Since in Eq. (88) we take into account all possible subsets, or partitions, we 
dropped there the use of the superscript. 
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quantifies the bipartite concurrence of the entangled part \<Pi2) of i.e., 
c^ 2 \^) = 0(^12), whereas and vanish identically for this bi-separable 
state. This selectivity with respect to correlations between specific subgroups 
of subsystems also hold for bi-separable mixed states - guaranteed by the 
construction (29) as convex roof. 

Another interesting quantity emerges for four-partite systems: Since the num- 
ber of subsystems is even, there is a term A^ 234 ^ = AP^ ® P^ ® P^ ® 
P-^ that does not contain any factor P+\ The corresponding concurrence 
vanishes for all states that do not contain proper four-partite correlations, 
i.e., bi-separable a nd tri-sep arable states. In particular, for a GHZ state, 
c^ 1234 ^(^ GHZ ) = 2^J2i^j K\j while c^ 23 ^ vanishes for H^-states - similarly 
like for tri-partite systems where H^-states bear only bipartite correlations . 
In the specific case of two-level systems, this particular choice of c^ 23 ^ mea- 
sures the potential of a given state for multi-particle teleportation [72] . 



3 Lower bounds 



In the case of bipartite two-level systems, we were able to evaluate the con- 
currence of arbitrary mixed states g exactly. This was possible because A (see 
Eq. (66)) was of rank one. For higher dimensional systems A typically is of 
higher rank, such that Eq. (76) exhibits two complications with respect to 
Eq. (72): additional non-linearities, and different matrices T a which, in gen- 
eral, cannot be diagonalised simultaneously. These two properties have so far 
prevented the derivation of an explicit solution of the optimisation problem 
formulated in Eq. (76). However, as we will show now, concurrence can be 
bounded from below, by some suitable approximations [30,31]. 



First, the Cauchy-Schwarz inequality [73] 



(89) 



allows to linearise Eq. (76). With x a 



N m 

c(g) > inf E y* 



=1 a=l 



yrpa e iip a yT 



yrpa e i<p a yl 



we conclude that 



> with E^« = 1 - 



(90) 



where we introduced some auxiliary real parameters y a . However, the different 
matrices T a still cause trouble finding the desired infimum. Here, the triangle 
inequality 

Ek|>|E^|> ( 91 ) 
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valid for arbitrary complex numbers z a , allows to circumvent this problem: 



For z a :=y a 



VT a e i(Pa V T 



, and y a > 0, one obtains 



N 



c(2)>mf£ 



i=i 



v(J2y a T a e^)v 



a=l 



(92) 



an expression for which the infimum is given analytically by Eq. (44). Our 
final expression for a lower bound of the concurrence therefore reads 



c(q) > maxjiSi - XJSj, °| ' 



with the singular values Sj of 

T = ]T Z a T a , Z a = y a e^, and ]T |Z a | 2 = 1 . 



(93) 



(94) 



The bound in Eq. (93) still depends on the choice of the Z a , what allows to 
tighten the estimate. Thus, one is left with an optimisation problem on an 2m- 
dimensional parameter space [30], where m is the number of matrices T a in 
Eq. (75). Note that the constraint J2 a \ Za\ 2 — 1 is by far simpler to implement 
than left unitarity, Eq. (31), since it is easily parametrised. Moreover, the 
dimension m of optimisation space is significantly reduced as compared to the 
dimension n\n\ of the original optimisation problem defined by Eqs. (29) and 
(31). 



3.1 Purely algebraic bounds 



The lower bound of Eq. (92) was obtained by the replacement of several ma- 
trices T a by a single suitably chosen T = J2 a Z a T a . So far, there is no clear 
prescription of how to chose the expansion coefficients Z a , which is partially 
due to the fact that the T a are determined only up to degeneracy. They are 
constructed with the help of the eigenvectors \x a ) of the projector A, what 
specifies an eigenspace, but does not distinguish states within these subspaces. 

One way to get rid of this ambiguity is to diagonalise A instead of A. For a 
typical state g, A will have no degenerate eigenvalues, and thus has uniquely 
defined eigenvectors T a . One can then provide m different lower bounds of 
concurrence directly calculating the singular values Sj in Eq. (93) for 

T = T a , a = l,...m. (95) 

In Section 3.3 we will see that one of these bounds alone may already yield a 
satisfactory approximation to the optimised lower bound. 
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3.2 Quasi-pure approximation 



By now, a large number of experiments focusses on the (controlled) prepa- 
ration and evolution of entangled states. Any degree of mixing - i.e. , the 
presence of classical correlations - decreases quantum correlations, and can 
lead to their complete destruction. Therefore - in particular in view of the 
many potential technical applications of non-classically correlated quantum 
states - pure entangled states are the experimentalist's desire: hence it is cru- 
cial to screen the investigated systems from the environment. 

In general, perfect screening is impossible, but experimental techniques are 
sufficiently advanced [5,8,74] to preserve entanglement over appreciable peri- 
ods of time. Yet, very little is known on the precise temporal evolution of such 
states even under weak but finite environment coupling, one of the principal 
obstacles being the lack of computable entanglement measures for arbitrary 
states. 

On the other hand, a general quantifier for entanglement, applicable to arbi- 
trary states, is not even needed in this context, since environmental influences 
can be assumed to be small under the given experimental conditions. Indeed 
the evolution of an initially pure into a mixed state occurs on a rather long 
time scale, and the experimentally interesting states are - though not exactly 
pure - at least quasi-pure, i.e. , they have one single eigenvalue /ii that is much 
larger than all the other ones. 

In order to provide some efficient means to deal with this type of problems, we 
now derive an analytic approximation of concurrence for quasi-pure states [31]. 
Indeed, we will find that this approximation also leads to a lower bound for 
arbitrary states. This quasi-pure approximation will allow for an efficient quan- 
titative treatment of non-classical correlations that arise in most present day 
experiments. 

The matrix A defined in Eq. (66) contains the matrix elements of A evaluated 
with the subnormalised (see Eq. (30)) eigenstates of g. Therefore, the elements 
of A are proportional to the eigenvalues of the considered state: 



Consequently, we can classify the elements of A according to their relative 
magnitude determined by the eigenvalues jij. This classification will serve 
as a basis for the approximate evaluation of concurrence in our subsequent 
treatment. 

The above proportionality leads to a natural order of the elements of A, in 
terms of powers of square roots of the real eigenvalues /ij of g, which we assume 




(96) 
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to be labeled in decreasing order, i.e. , [i\ >> [ii > . . . > fi n . Hence, if we con- 
sider terms proportional to either one of the fij, with j > 1, as perturbations 
of the dominant term A\\ ~ /if, we obtain the following classification: 

- the element A\\ is lowest order, 

- all elements with one index j > 1, i.e. Aj\, A\j, A{\ and A\{, are first order, 
and 

- elements with two indices j, k > 1, alike A 1 ^. or Ajl, are second order. 

However, this classification is not yet a sufficient basis for our approximation. 
In fact, the element A\\ is lowest order, though still it could vanish. This is the 
case if and only if the sub normalised eigenstate \ipi) to the largest eigenvalue 
/ii of q is separable, since A\\ is the square of the concurrence of j^i) - see 
Eqs. (62) and (66). Therefore, as an additional requirement to quasi-purity 
we have to impose that j^i) is entangled. Since the desired approximation 
is supposed to be applied to states that occur in the experiments mentioned 
above, this is not too stringent a restriction - if an ideal experiment without any 
environment coupling led to a pure state with non-negligible entanglement, it 
is reasonable to assume that the eigenstate j^i) associated with the largest 
eigenvalue /ii of g is not separable either. 

We wish to approximate A by a matrix product 

Af k ~ T* m T jk , (97) 

with a complex symmetric matrix T e C nxn . Such a replacement allows for 
an analytic solution, since the sum over a in Eq. (76) reduces to a single term, 
and the analytic expression Eq. (39) for the infimum, derived in Section 2.1, 
can be employed. 



Specifically for the lowest order term A\\, Eq. (97) yields 7n = yA\\, up to an 
arbitrary phase which can be dropped. Subsequently, evaluation of Eq. (97) 

for the first order elements leads to 7^ = Ajl/\jA{\, employing also the 
expression for 7n. Finally, we still have the freedom to fix Tj k , for j,k ^ 1. 

For this purpose we use Eq. (97) for A^l, what leads to T jk = A 1 ^/ \J A\\, such 
that all matrix elements of T are given by 

A 11 

T jk = -7== . (98) 
An 

With this choice of T, Eq. (97) is exact at lowest and first order, and - in 
addition - the second order elements Aj k are taken into account correctly. Note 
that, using only one single matrix T, it is not possible to describe accurately 
all second order elements, such as A*\. All terms of third and fourth order are 
dropped in the present approximation. 
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Since we were assuming j^i) to be entangled, i.e. , A\\ finite by virtue of 
Eq. (62), T is well defined. Approximating A in terms of this matrix T, 
Eq. (76) can be approximated as 



c(g)^c qp (g) = mfY,\[vrV T } u \ , (99) 



following the same steps as in the derivation of Eq. (72). As formulated in 
Eq. (44), the infimum can be expressed in terms of the decreasingly ordered 
singular values <Sj of T, 



c, 



qp 



(g)=max{s 1 -^2S i , . (100) 



i>l 



A priori, it is not clear up to which degree of mixing the quasi-pure approxi- 
mation provides reliable estimates. The fact that not all second order elements 
are taken into account may somewhat reduce our confidence in a wide range 
of applicability for this approximation. However, we will see in Sections 4.1, 
4.2.3, and 4.2.4 that the above ansatz provides very good results for many 
states, even for those with a substantial degree of mixing. 

Indeed, the quasi-pure approximation is not only an approximation, but even 
a lower bound of concurrence: Given the decomposition A l j% = Z)a(^]m)*^fe' 
T is easily expressed as 

T = Y,z*T a , with z a = , (7ri) . (101) 

v/E^riT 

Thus, T is indeed a valid symmetric matrix to provide a lower bound as 
formulated in Eq. (92). 



3.3 Lower bounds of states with positive partial transpose 



In the preceding section we collected a set of operational tools to access the 
non-classical correlations of arbitrary mixed quantum states of finite dimension 
- characterised by their concurrence. Our formal treatment spans the entire 
range from an approximation-free description for numerical implementation, 
over lower bounds that can be tightened numerically, to an easily tractable, 
purely algebraic estimate of the degree of entanglement of a quasi-pure quan- 
tum state which is typically dealt with in experiments. However, our approach 
is largely based on physical intuition, and so far we cannot come up with math- 
ematically precise error bounds. Often, however, the latter are available in full 
mathematical rigour only under rather restrictive assumptions - while we are 
seeking for robust quantities which, beyond formal consistency, can cope with 
requirements which stem from real-world experiments. 
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Fig. 4. Numerically optimised lower bound (dashed line, Eqs. (93) and (94)) of the 
concurrence of the family of bipartite spin-1 states g a [75] as defined in Eq. (102), 
together with the algebraic bounds (solid lines) obtained from Eq. (93) using (95). 
Both the optimised bound and the largest algebraic bound are positive, such that 
the state is detected as entangled in the entire parameter range a = [0, 1]. All other 
algebraic lower bounds are negative. Note that the scale of concurrence (y-axis) is 
different for positive and negative bounds. 

Now we will test our toolbox under realistic conditions, and we start out with 
the detection of nonseparable states with positive partial transpose (ppt). In 
the next chapter, we will then use our approach to monitor the time evolution 
of entanglement, under various scenarios. 



3. 3. 1 Some exemplary ppt states 

One of the main requirements imposed on an entanglement measure is that it 
be able to distinguish entangled states from separable ones. Whereas a large 
class of entangled states is detected by the ppt criterion denned in Eq. (12), 
no operational criterion is known so far that can detect all states with positive 
partial transpose. In general it tends to be rather demanding to decide whether 
such states are entangled or not. Therefore - albeit our bound, Eq. (92), is 
capable of more than just checking separability - we use it, as a first test of 
its pertinence, as a separability criterion for some families of entangled states 
with positive partial transpose [39,40,75]. 

The first class of states describes a bipartite spin-1 system. The state g a acting 
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Fig. 5. Algebraic lower bounds of the concurrence (see Eqs. (93) and (95), solid 
lines) of a family of 4 x 2 states [75] denned by Eq. (103), plotted as a function of 
the parameter a. Although none of these bounds is positive, the optimised bound 
(dashed line), Eq. (93), is positive. Thus the state is detected as entangled in the 
entire parameter range of a. 

on C 3 <g> C 3 is given [39,75], for a G [0, 1], as 



Qa = 



1 + 8a 
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(102) 



and has a positive partial transpose as denned in Eq. (12), in the entire range 
of a. The algebraic lower bounds denned by Eqs. (93) and (95) are plotted in 
Fig. 4 as solid lines. One of them is positive for all values of the parameter 
a, and the non-separability of g a is therefore detected by a purely algebraic 
criterion. All other algebraic bounds are negative, and therefore do not provide 
any information on their own. 

The dashed line in Fig. 4 shows the lower bound that is numerically optimised 
over the Z a in Eq. (92), using a downhill simplex method [76]. It is significantly 
larger than the positive algebraic bound, and shows a qualitatively different 
behaviour for large a, where its first derivative is finite, while that of the largest 
algebraic bound vanishes for a = 1. 

A second class of states g a [75] acts on C 4 Cg>C 2 , and once again has a positive 
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a 

Fig. 6. Algebraic lower bounds (see Eqs. (93) and (95), solid lines) of a family 
of 3 x 3 states [40] defined in Eq. (104), as a function of the parameter a. For 
a < 1/2, g a is separable, and for a > 1/2 it is entangled [40]. It has positive partial 
transpose for a < 3/2, and non-positive partial transpose for a > 3/2. The dashed 
line shows the numerically optimised lower bound, Eq. (93). The state g a is detected 
as entangled by the algebraic bound, exactly in that parameter range where it has 
non-positive partial transpose. It is detected as entangled by the optimised bound 
in approximately half the parameter range with positive partial transpose, what, 
however, might be due to a failure of our numerical optimisation routine. 

partial transpose, Eq. (12), for a G [0, 1]: 



Qa 



a0000a00 
OaOOOOaO 
OOaOOOOa 
1 OOOaOOOO 
l + 7 a 0000/3007 
aOOOOaOO 
OaOOOOaO 
00a0700/3 



with 



1 + a 

1 2 ' 

7 = 7i 



(103) 



Figure 5 shows the algebraic lower bounds obtained as discussed in Section 3.1, 
which are all negative. Thus, none of them detects g a as entangled. Though, 
due to the degeneracy of this particular state, there is a degeneracy in the 
eigenvalues of A, such that the matrices T a , and, consequently, also the al- 
gebraic lower bounds, are not uniquely determined. Neither did we find any 
matrices T a in the degenerate subspaces that provide positive lower bounds. 
Yet, the numerically optimised lower bound - also shown in Fig. 5 - is positive 
in the entire parameter range. Hence, also this state is detected as entangled 
by our lower bound, Eq. (93). 
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A third class of states g a [40], acting on C 3 ®C 3 , is defined for a G [-5/2, 5/2], 
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with 



P ± = -±a 



(104) 



Since replacing a by —a is equivalent to exchanging the subsystems, we will 
discuss this state only for a G [0, 5/2]. The state g a has a non-positive partial 
transpose for a G [3/2,5/2], is entangled with positive partial transpose for 
a G [1/2, 3/2] and is separable for a G [0, 1/2] [40]. As depicted in Fig. 6, g a is 
detected as entangled in its domain of negative partial transpose already by 
the best algebraic lower bound. In the regime where g a has positive partial 
transpose all algebraic bounds are negative, such that the optimised lower 
bound is required for distinguishing g a from separable states. However, even 
the optimised bound does not detect g a in the entire interval a G [1/2, 3/2]. For 
a < 1.02, the lower bound seems to fail as a sufficient separability criterion. 
At present, we have no conclusive evidence from our numerical optimisation- 
routine to decide whether the bound itself is not good enough, or whether the 
numerically found maximum is not the global one. 

The above exemplary ppt states show that our lower bound, Eq. (93), is capa- 
ble to detect families of entangled states which are not recognised by the ppt 
criterion. For some states it is even not necessary to evaluate the optimised 
bound, since already one of the algebraic bounds, introduced in Section 3.1, is 
positive. Moreover, also the quasi-pure approximation is positive for some ppt 
states [31]. However, there are also states with only negative algebraic bounds 
and negative quasi-pure approximation, though positive optimised bound. Our 
last example above showed a case of entangled states that we have so far been 
unable to detect for a small subset of parameters, though it remains hitherto 
undecided whether this is a failure of our numerical optimisation routine or 
of our lower bound, Eq. (93), itself. It is obvious from the different behaviour 
of optimised and algebraic bounds at the border line between ppt and non-ppt 
regions in Fig. 6, that some more profound algebraic signatures remain to be 
uncovered. 



4 Dynamics of entanglement under environment coupling 



Arguably the central motivation for deriving efficiently evaluable measures of 
the entanglement of mixed states is the ubiquity of the latter in any real physi- 
cal setting. If we consider entanglement as the central resource of most types of 
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quantum information processing, then the experimentally most relevant ques- 
tion is that of the lifetime of entanglement under the environment-induced 
mixing. This is the subject of the present, concluding section of this review. 

In a first subsection, we will test the mixed state entanglement estimates 
derived in the previous sections, such as to demonstrate their versatility to 
describe generic time evolutions under environment coupling. Here, the time 
evolution both of the system and of the bath will be generated by random 
Hamiltonians, without any specific physical realization in mind. 

Subsequently, we will specialise to particular, experimentally relevant (since 
realized) cases, and specifically focus on the scaling properties of the typical 
time scales which determine the time evolution of bipartite and multipartite 
entanglement. 

4-1 Random time evolution of higher- dimensional bipartite systems 

Let us first have a closer look at the performance of the various entangle- 
ment estimates derived above, for a generic time evolution under environment 
coupling. For that purpose, we consider a bipartite system and a third sys- 
tem serving as environment. The bipartite system is initially prepared in a 
maximally entangled pure state 

l*o> = 4=X>> ' (1 ° 5) 

yd i=1 

i.e. , it is not entangled with the environment. We then evolve the total system 
under a unitary dynamics generated by the Hamiltonian 

H = H sc + H s <8> l c + l s <g> H e , (106) 

with a randomly chosen hermitian matrix H se acting on system and environ- 
ment, while a second randomly chosen hermitian matrix H s only acts on the 
(bipartite) system, but not on the environment. All elements of H sc and H s are 
determined independently under the constraint of hermiticity. Any real entry 
is obtained as sin(r), with a random integer r, where the a priori probability 
is the same for any integer < r < 10 15 [77]. Furthermore, we neglect the free 
evolution of the environment by setting H e = l e . 

Tracing out the environment after a finte interaction time t leads to a mixed 
state of the bipartite system. H s describes the interaction between the bipartite 
system's components and does not induce any mixing, though it can change 
the degree of entanglement. We further scale the Hamiltonians H se and H s by 
the real coupling constants a sc and a s , which fix the relative time scales of 
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Fig. 7. Top panel: Upper bound (Eq. (79), solid line) and numerically optimised 
lower bound (Eqs. (93) and (94), dashed line) of concurrence of a 3 x 3 state, vs. 
scaled time a s t. The system is initially prepared in a maximally entangled state of 
type (105), with concurrence c = 2/y/3 (see Eq. (61)). A random time evolution 
according to (106), with a sc /a s = 10~ 2 , leads to a finite degree of mixing measured 
by the system state's von Neumann entropy (dash-dotted line). The largest eigen- 
value of g evolves along the dash-double-dotted line. In the lower panel, the solid 
and the dashed line represent the difference between the optimized lower bound of 
concurrence, Eqs. (93) and (94), and the quasi-pure approximation, Eq. (99), or the 
best algebraic lower bound, Eqs. (93) and (95), respectively. As also visible from 
the upper panel, all our lower bounds provide excellent estimates of the actual value 
of concurrence over the entire evolution period, almost undistinguishable from the 
upper bound, Eq. (79) (full line in the upper panel). 

unitary and incoherent system dynamics. a se determines the strength of the 
system-environment interaction, and therefore the mixing rate of the system, 
whereas a s specifies the time-scale of the unitary system evolution, causing a 
reversible decrease and increase of bipartite entanglement therein. 

With these premises, we can now monitor the time evolution of the concur- 
rence under the above random dynamics, as a function of the scaled time a s t, 
for different values of the ratio a sc /a s . The degree of mixing of g will be char- 
acterised by the von Neumann entropy S = —Tig In g of the system state. In 
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addition, we also follow the time evolution of the largest eigenvalue of g, since 
this was the central pillar in our derivation of the quasi-pure approximation 
(Section 3.2). 

The upper panels in Figs. 7 and 8 show the time evolution of upper (Eq. (79), 
solid line) and lower (Eqs. (93) and (94), dashed line) bounds of concurrence, 
together with concurrence in quasi-pure approximation (Eq. (99), dotted line), 
for two different random time evolutions of the same maximally entangled 
3x3 initial state. The initial value of c(g) follows immediately from Eq. (61). 
The interaction Hamiltonians H sc and H s are the same in both figures, but the 
coupling constant a se in Fig. 8 is by a factor five larger as compared to its value 
a sc /a s = 10~ 2 in Fig. 7. Consequently, mixing increases slower in Fig. 7 than 
in Fig. 8. The degree of mixing is characterised by both the largest eigenvalue 
of g, depicted as a dash-double-dotted line, and by its von Neumann entropy 
(dash-dotted). 

For both values of a se /a s , the quasi-pure approximation and the optimized 
lower bound remain almost indistinguishable, over the entire system-environment 
interaction time. The difference of both quantities is plotted in the lower panels 
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Fig. 9. Same as Fig. 7, for a 4 x 4 system initially prepared in the maximally 
entangled state, Eq. (105), and a se /a s = 10~ 2 . 



of the figures, together with the difference between the numerically optimised 
lower bound, Eqs. (93) and (94), and the best algebraic bound, Eqs. (93) 
and (95). Over almost the entire time interval displayed in Fig. 7, and for 
short times in Fig. 8, where mixing is not too large, the relative error of both 
approximations is about two orders of magnitude smaller than the actual value 
of concurrence itself (note the different scales on the vertical axes). Only states 
with rather small concurrence and a large degree of mixing exhibit a significant 
difference between the optimized lower bound and the best algebraic bound 
or the quasi-pure approximation. 

Fig. 9 repeats the scenario of Fig. 7, for a 4 x 4 system, with environment cou- 
pling constant a sc /a s = 10~ 2 . The degree of mixing remains moderate during 
the entire time evolution, as indicated by the von Neumann entropy and the 
largest eigenvalue of g. Once again, in this higher-dimensional system, upper, 
Eq. (79), and optimised lower bound, Eq. (93) and (94), provide an excellent 
estimation of the actual value of concurrence. Equally so, also the best al- 
gebraic lower bound, Eqs. (93) and (95), and the quasi-pure approximation 
perform very well, over the entire time interval. 
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Fig. 10. Same as in Figs. 7 and 9, for a 2 x 4 system initially prepared in the 
maximally entangled state, Eq. (105), and a sc /a s = 10~ 2 . 

Qualitatively the same observations are made for the time evolution of an 
initially maximally entangled state, Eq. (105), of a 2 x 4 system, monitored in 
Figs. 10 and 11, for different system-environment coupling strengths a sc /a s = 
10~ 2 and 5 x 10~ 2 , respectively. Only during a rather short time interval, when 
concurrence exhibits a dip at a s t ~ 3 x 10~ 4 in Fig. 11 does the discrepancy 
between optimal lower bound and its estimates increase to about one percent 
of the actual value of concurrence, which, however, rather highlights the ex- 
cellent reliability of our estimates for all practical purposes. In particular, this 
holds true for the quasi-pure approximation, albeit the strongly mixed states 
occuring as time evolves in Fig. 11 hardly satisfy the basic assumptions made 
for its derivation (see Section 3.2). 

4-2 Realistic scenarios of entanglement dynamics 

Recently, many groups were able to prepare entangled states in a variety of 
physical systems and experimental setups, demonstrating an impressive ability 
to manipulate and detect them efficiently [7,8,9,10,11,12,13,14,15,16,17]. Par- 
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Fig. 11. Same as Fig. 10, but for a larger system-environment coupling 
a sc /a s = 5 x 10~ 2 . 

ticular effort has been devoted to entangle multiple subsystems, not only to 
investigate the possibility of fulfilling the fundamentally important scalability 
requirements for quantum computation, but also to understand how far one 
can push a quantum system towards the macroscopic limit and still observe 
entanglement - an intrinsic quantum feature with no classical counterpart. 

However, a major obstacle for the controlled entanglement of more and more 
subsystems remains with the incapacity of achieving perfect screening of the 
system from the environment. After some time, the unavoidable residual in- 
teraction with the reservoir induces mixing of the system state, and thus the 
emergence of classical correlations, at the expense of quantum entanglement. 
Hence, we face the highly relevant task of understanding the sources of entan- 
glement decay, what implies the identification of the associated time scales. 

Although a general solution to this problem, for arbitrary system dynamics 
and decoherence mechanisms is still out of reach, our technical machinery 
developed in the previous sections allows to treat arguably all situations en- 
countered in typical state of the art experiments, as we shall outline in the 
sequel. We start with a short recollection of environment models adapted for 
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decoherence processes in a typical quantum optical context. 



4-2.1 Environment models 

We already anticipated in Section 4.1 that an open system can be described as 
a system S which interacts with an environment E, such that the dynamics of 
the total system S + E is unitary, and governed by the von- Neumann equation 

^ = 4 [*.*-!. <«") 

with the total Hamiltonian of the form (106). However, we are only interested 
in the properties of the system S itself. An equation for the evolution of the 
(reduced) system density matrix g is obtained by performing a trace over the 
environmental degrees of freedom: 

g = *J- (108) 

In Section 4.1 we considered S as a bipartite system and explicitely performed 
the total evolution given by H, obtaining the reduced density matrix after 
tracing over environmental degrees of freedom. However, an equation only for 
the density matrix of the system, after an interaction time t, can be obtained. 
Under the assumption of complete positivity and Markovian dynamics, it can 
be written in the Lindblad form [78,79,80] 

^ = — % - [H s , q] + Cq= ~ [H s , g] + y ( 2 d i sA - d\ dig- gd\ d^j , 

(109) 

where operators describe the system-environment coupling, with strength 
ry The use of the Markovian approximation is well justified in a large vari- 
ety of quantum optical experiments where entanglement has been produced, 
although one should mention that non-Markovian effects can be important in 
the description of some condensed-matter systems [81]. 

Different situations may arise when a system is coupled to the environment: 
energy can be exchanged and dissipation can take place, noise can be added to 
the system, or elastic processes can introduce loss of phase coherence without 
energy transfer. All these processes can be described in terms of the above 
master equation by a suitable choice of the operators d iy which can be written, 
in the case of two-level systems, in terms of the Pauli matrices. For a two-level 
system interacting with a thermal bath, for example, the non-unitary part of 
the master equation reads [82] 

T(n + 1) , n . Tn , n 
Lq = (2 a_ga + — a + o-_Q — ga + a_)-\ (2 a + ga_ — a_a + g — ga_a + ) 

(110) 
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In this equation, the first and the second term on the right hand side describe, 
respectively, decay and excitation processes, with rates which depend on the 
temperature, here parametrised by n, the average thermal excitation of the 
reservoir. In the limit of vanishing temperature, n — 0, only the spontaneous 
decay term survives, leading to a purely dissipative process, 



which drives the system, asymptotically, to its ground state. 

Another important limiting case of equation (110), which describes noisy dy- 
namics, is obtained for infinite temperature, where n — > oo, and, simultane- 
ously, r — > 0, so that Tn = F remains constant: 

f f 

Cg — — (2 a^ga + — a + a^g — ga + a ) + — (2 a + ga^ — a^a + g — ga-a + ) . 



In this case, decay and excitation occur at exactly the same rate, and the 
noise induced by the transitions between the two levels brings the system to 
a stationary, maximally mixed state. 

A purely dephasing reservoir is obtained by choosing di = d = a + a_ in 
Eq. (109) leading to the master equation 



In this case, there are no changes in the populations of the ground and excited 
states, and energy is conserved. Only the off-diagonal elements of the reduced 
density matrix decay, leading to a loss of phase coherence. 

In the context of entanglement, we have to describe how the environment 
acts on the composite system. In Sections 4.2.2 and 4.2.3 we will assume that 
each subsystem interacts independently with the environment, a well justi- 
fied assumption whenever the particles composing your system are sufficiently 
separated from each other and, therefore, no collective environment effects 
must be taken into account. Hence, each particle is subjected to a dynam- 
ics described by its own Lindblad operator, supposed to be of the same form 
for all the components of the system. Note that this assumption of mutu- 
ally independent environment couplings immediately implies that all initially 
entangled states will asymptotically evolve into a separable state. We shall 
therefore focus on the precise timescales of this decay process. 

In the final Section 4.2.4, we will consider an entanglement scheme [83] which 
has been experimentally implemented using ion traps [9], and in which deco- 
herence acts indirectly through the coupling of the collective center-of-mass 



dg _ T 
~dt ~ 2 



(2cr_£>(7 + — a + a-g — ga + a^) , 



(111) 



(112) 




(113) 
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motion of the ions with the environment. Here, again, environment interaction 
has a detrimental effect, which we shall study in more detail focusing on the 
scaling of multipartite entanglement with the system size (i.e., the number of 
entangled particles). 



4-2.2 Entanglement dynamics of bipartite two-level systems 

For bipartite two level systems, one can use the exact expression, Eq. (39), for 
concurrence, to follow the time evolution of entanglement. While much simpler 
than the multipartite case, the bipartite situation is useful for developing some 
intuition about the entanglement decay under different environment dynamics. 

To start with, we consider an initial Bell state |^ ± ) = (|01) ± |10)) /y/2. 
The solution of the master equation for the different environment couplings 
discussed previously is straightforward, and the time evolution of concurrence, 
shown in Fig. 12, is given by 

f e -4ft . i ) 

c(t) = max I — + e~ 2rt - - , , (114) 

for the infinite temperature case, and by 

c(t) = e~ rt , (115) 

for dephasing, Eq. (113), and zero temperature, Eq. (Ill), environments. The 
situation is the same for the Bell state |$ ± ) = (|00) ± |11)) /V2 (right panel of 
Fig. 12), apart from the zero temperature case where the concurrence decays 
as c(t) = e~ 2Tt . This accelerated (by a factor of two) decay of concurrence for 
the states as compared to the l^/ 1 * 1 ) states, under the influence of zero 
temperature environment, can be understood from the time scales involved in 
the corresponding solution for the density matrix: while for \^) each term 
1 1 ) and 1 10) corresponds to a single particle decay, leading to a time scale 
e~ rt , we have the term 1 11) in l^), such that both particles can undergo an 
environment induced transition to the ground state, thus introducing a faster, 
e~ 2rt , decay. 6 General solutions for the finite temperature case can be given 
explicitely [84], though are rather cumbersome and will not be presented here. 
However, one may expand these solutions and obtains, at first order in t, 

c(^ ± , t) ~ 1 - (2n + 1 + 2^n(n + 1)) Tt, (116) 



6 Although intuitive and reasonable in this case, one must be carefull while using 
such kind of arguments to deduce the actual behaviour of concurrence from the 
time scales which appear in the coherences of g. In general, there is no simple and 
obvious relation between these quantities. 
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Fig. 12. Time dependence of the concurrence of a bipartite two-level system initially 
prepared in Bell states |^ ,± ) (left panel) or |$ ± ) (right panel), under coupling to a 
zero temperature (dotted lines), dephasing (long dashed lines), infinite (solid lines), 
and finite temperature (n = 0.1, dashed lines) environment (time evolution gener- 
ated by Eqs. (Ill), (113), (112), and (110), respectively. In the case of states, 
the zero temperature and dephasing solutions coincide. In contrast, the zero tem- 
perature environment induces a twice as fast decay as compared to the dephasing 
environment, for the \^} state. For all positive temperatures these initial states 
evolve into separable states in a finite time. Only zero temperature and dephasing 
environments induce separability only in the limit t — > oo. 



and 

c($ ± ,t) ~ l-2(2n + l)rt, (117) 

for and \^) states, respectively. These expressions not only reproduce 
the leading order terms of the zero and infinite temperature solutions in the 
appropriate limits, but also show the influence of the temperature on the short 
time behaviour of concurrence decay. 

It is equally simple to obtain the long-time asymptotics for the singular values 
of the matrix r defined in Section 2, and, consequently, for the quantity that 
enters expression (39) for concurrence. It reads, both for and states, 



•°°V U 7 (2« + i) 2 



This quantity is non-positive for all values of n and vanishes only in the case 
of zero temperature. This means that the above initial states evolve into sep- 
arable states within a finite time, for any T > 0. In contrast, for zero temper- 
ature as well as for dephasing environments does c(t) only vanish in the limit 
t -> oo. 7 



Note that the above is not true for general states. In fact, some initial mixed states 
may reach separability on finite time scales also in the zero temperature case [85]. 
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4-2.3 Entanglement dynamics of multipartite two-level systems 

Let us now generalise our above investigation for larger arrays of two-level 
systems, with variable size N. These are objects often encountered, e.g., in 
ion trap quantum computing schemes, and the robustness of their entangle- 
ment properties with increasing size is at the very heart of the fundamental 
scalability requirement on which hinges any future technological application. 

We have already prepared a comfortable tool for such generalisation - the 
multipartite concurrence defined in Eq. (88), with the particularly useful 
property 



c N (\<P)® \<p)) = c N - 1 (\$)) , (119) 
for \$) <g) \tp) e Hx <g> . . . <8> H N , and \tp) e H h i E {1, . . . , N} . (120) 

This allows for the direct comparison of the entanglement properties of pure 
and mixed states of multipartite qubit arrays of increasing size [32], where 
the mixed state concurrence is once again evaluated through the convex roof 
Eq. (67), with the appropriate definition of the operator A in Eq. (87). 

Since the reliability of the lower bounds derived in Section 3 was demonstrated 
extensively in Section 4.1 above, we will actually use the quasi-pure approx- 
imation Eq. (99) in the following, minimising the numerical effort very con- 
siderably, notably for large system sizes N (notwithstanding, some additional 
tests were performed at randomly chosen instances, to convince ourselves of 
the quality of the quasi-pure approximation with respect to the optimal upper 
and lower bounds, Eqs. (79) and (93), respectively. 

We shall consider two types of initial states, somewhat similar to the singlet 
and triplet states in Section 4.2.2 above: the GHZ state 

|^v)ghz = -^(|00...0) + |ll...l)) , (121) 

and the W state, 

\V N )w = A= (|00 . . . 01) + 1 00 . . . 10) + . . . + 1 10 . . . 00)) . (122) 
v N 

These were recently produced [9,12] in the lab, for N — 3, 4, and it is now 
within experimental reach to monitor the time dependence of their degree of 
entanglement, by means of quantum state tomography [17]. Figure 13 shows 
the decay of the concurrence for tri-partite GHZ and W states under the 
influence of zero temperature (circles), dephasing (triangles) and infinite tem- 
perature (squares) environments as a function of time. The symbols represent 
the result of the quasi-pure approximation, while the lines indicate the best 
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Fig. 13. Time dependence of the multipartite concurrence c/v(i), ^q. (88), for N = 3 
particle systems initially prepared in a GHZ (dashed lines) or W (solid lines) state, 
under the influence of a zero temperature (circles), dephasing (triangles), and infinite 
temperature (squares) environments. In all cases, the numerical results are very well 
fitted by an exponential decay (solid and dashed lines). 

fit to a mono-exponential decay cjv(i) = Ae~ l1 + B. As in the bipartite case, 
the zero temperature and dephasing environment lead to separability only in 
the limit t — > oo, for these initial states. 8 

To assess the scaling properties of multipartite entanglement for variable sys- 
tem size, we now compute the time evolution of concurrence, precisely as in 
Fig. 13, for different values of N, and extract the decay rates 7. Fig. 14 shows 
the scaling of the entanglement decay rates 7 with N, under the above, distinct 
decoherence mechanisms. 

We see that the decay rate of the concurrence of GHZ states (top panel of the 
figure) increases linearly with N, except for the small- N behaviour of 7 for 
the zero temperature environment. Indeed, in the special case of dephasing 
environment, the density matrix is always a mixture of two pure states, and, 
hence, of rank two, and can be treated analytically. The observed behaviour 
comforts our intuition - which suggests that the larger the system, the easier it 
is for detrimental environment effects to manifest. In addition, this fragility of 
GHZ states was also observed in [86] by an analysis based on their separability 
and distillability properties [87,88] rather than a dynamical approach. 

Remarkably, the situation changes quite drastically for the W states (bottom 
plot of Fig. 14). In this case, only the infinite temperature environment gives 
rise to an almost linear increase of 7 with N, slightly faster than for the GHZ 

8 Note that, as in the bipartite case, this is not true for general states. The state 
(| 000) + 1 01 1) + 1 101) + 1 110) ) /2, for example, reaches separability, according to 
quasi-pure calculations, on a finite time scale also in the zero temperature case. 
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states. In contrast, for dephasing and zero temperature reservoirs, the decay 
of the concurrence is independent of N. Moreover, the zero temperature case 
also allows for an analytic solution (as above, the rank of the state reduces to 
two) for all N, leading to c N (t) ~ e~ rt . Consequently, for these environments, 
W states clearly outperform GHZ states in terms of the robustness of their 
multipartite entanglement properties. 9 

One might be tempted to attribute this to the smaller initial concurrence of 
W as compared to GHZ states (see Fig. 13), though the ratio 



with a maximum at N = 5, approaches unity for large N (with a value 1.07 
for N = 7, the maximal size considered in Fig. 14). 

4-2.4 Experimental entanglement production: the M0lmer-S0rensen scheme 

In the previous sections we assumed that a perfect pure entangled state is 
available at the beginning, and monitored the reservoir induced entanglement 
decay in the course of time. Here, we will consider a situation where the en- 
vironment acts simultaneously to a unitary evolution which is intended to 
prepare the entangled state. Specifically, we scrutinize a scalable scheme pro- 
posed in [83] to produce GHZ-like states (see Eq. (121)), 



experimentally implemented at NIST [9] to produce controlled entanglement 
of two and four trapped ions. 

To illustrate this preparation scheme - which is valid for iV ions - let us 
consider the case N = 2 with two ions confined by a harmonic potential, si- 
multaneously illuminated by two electromagnetic fields. The energy levels are 
depicted in Fig. (15), where uq is the frequency of the electronic transition, 
and v the oscillation frequency of a given collective motional mode of the 
particles in the trap. The fields oscillate with c^o + v — A and ujq — v + A in 
such a way that the two photon process that drives the 1 00) <-> |11) transition 
is resonant and a superposition of these states can be produced. At the end 
of the preparation scheme, all one-photon processes, which excite motional 
states, interfere destructively with each other, such that vibrational and in- 
ternal degrees of freedom are uncorrelated. However, during the illumination 
these correlations do exist, and motional decoherence can decrease the success 

9 Cluster states also present a kind of robust behaviour in terms of separability and 
distillability criteria in the case of a depolarising channel [86]. 




tf JV ) = -^(|00...0)+e*"|ll...l)) , 



(123) 
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Fig. 14. Concurrence decay rates 7 (in units of the reservoir rate T) for GHZ (top) 
and W (bottom) states, as a function of the system size (i.e., the particle number) 
N. The different environment models are represented by circles (zero temperature), 
squares (infinite temperature), and triangles (dephasing). The connecting lines are 
drawn to guide the eye. Whilst for GHZ states the decay rates increase roughly lin- 
early with N (for sufficiently large N), independently of the specific environment, 
the W states exhibit increasing decay with system size only for the infinite temper- 
ature environment. Remarkably, the decay rate of the W states is size-independent 
for dephasing and zero temperature environments! 

probability of the preparation process. The decay rate of the electronic states 
is negligible during the time scale of the state preparation, and therefore one of 
the main sources of errors in the system arises indirectly through the coupling 
with these decohering motional modes. 

The heating of the ions, related to fluctuating fields in the trap electrodes, 
leads to a thermal motion with steadily increasing temperature and can be well 
described by the infinite temperature reservoir discussed previously, Eq. (112), 
with measured heating rates Y/u ranging from 10~ 4 to 1CT 3 [89,90]. 

Figure 16 shows the time evolution of the multipartite concurrence cat, Eq. (88), 
(evaluated through its quasi-pure approximation, Eq. (99)) under this scheme, 
for four ions and heating rates Yjv = (bold solid line), 1 x 10~ 4 (solid line), 
2 x 10" 4 (dashed line), 3 x 1CT 4 (dot-dashed line) and 4 x 1CT 4 (dotted line). 
The system starts with all ions in the electronic ground states and, during 
the evolution, entanglement builds up until reaching its maximum value at 
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Fig. 15. Illustration of the M0lmer-S0rensen scheme [83] for the creation of a max- 
imally entangled state (|00) + i\ll))/y/2. Two ions are illuminated simultaneously 
by two electro-magnetic fields. One (dashed arrows) is red detuned with respect 
to the blue sideband, with detuning A. The second pump field (dotted arrows) is 
blue detuned with respect to the red sideband. Both sideband transitions are not 
driven resonantly, however the two photon process 1 00) <-> |11) is resonant, since the 
absolute values of the detunings of both transitions coincide. Thus, starting with 
the initial state |00), one can create a coherent superposition of the states 1 00) and 
|11). 

vtp re p — 7r/ ~ 1100 (with f] the Lamb-Dicke parameter, the ratio between 
the width of the vibrational ground state and the wavelength of the driving ra- 
diation, and Q the single ion resonant Rabi frequency describing the strength 
of ion- field coupling). After this preparation time the lasers should be turned 
off and, for zero environment coupling, concurrence achieves its ideal value 
corresponding to the GHZ-like state, while it shows appreciable decrease even 
for small heating rates. 



As in the case of the purely environment induced dynamics, we are interested 
in the scaling of the achieved degree of entanglement with the system size N. 
In the ideal case, the maximum value of the multipartite concurrence, Eq. (88), 
increases with N as = 2 1 ~ JV / 2 y / (2 iV — 2)/2, saturating at \[2 for large N 
as depicted by the bold solid line and circles in Fig. 17. With the addition 
of motional heating, this growth with N is not monotonic anymore and we 
observe, for example, that from a certain value of the coupling strength the 
maximum concurrence for six ions gets smaller than the one for five. Hence, 
the detrimental effect of (indirect) motional decoherence for entanglement 
generation is enhanced with increasing system size. This is further illustrated 
in Fig. 18, where we plot the entanglement loss A c = c% HZ — c™ ax , i.e., the 
difference between the ideal value and the actual optimal value achieved at a 
finite decoherence level. For all values of the ion heating rate, the entanglement 
loss increases with the system size. 10 



Although one might be tempted to associate this behaviour with the fragility of 
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Fig. 16. Quasi-pure calculation (see Section 3.2) of the time evolution of the mul- 
tipartite concurrence for four ions. In the case of vanishing environment cou- 
pling (bold solid line), concurrence reaches its ideal value at the preparation time 
z^prep = 7r /( r /^) — 1100, with r\ = 0.05 and O, = 0.057. As we increase the cou- 
pling strength T/u, the maximally achievable entanglement decreases, as shown for 
T/v = 1 x 10~ 4 (solid line), 2 x 10" 4 (dashed line), 3 x 10~ 4 (dot-dashed line) and 

4 x 10~ 4 (dotted line). 

5 Conclusion 

In the present contribution, we have performed a program reaching from the 
formal definition to the time dependent monitoring of concurrence of arbitrary, 
finite dimensional, multipartite quantum states under incoherent environment 
coupling. The general applicability of the various lower bounds of concurrence 
which have been derived stems from the algebraic structure of a suitable re- 
formulation of its original definition, with direct generalizations for higher 
dimensional and/or multipartite quantum states. 

We have checked the reliability and tightness of our various approximate lower 
bounds of concurrence under very different physical scenarios, from rather 
untypical entangled states with positive partial transpose to experimental 
schemes of controlled environment production. Whilst these results generally 
suggest a rather excellent performance of the here derived, efficiently com- 
putable entanglement measures, also quite a few questions remain to be an- 
swered on the mathematical side: Which of our lower bounds of concurrence 
are actually able to detect all entangled states? Are there analytical error 

GHZ states with increasing number of particles as discussed in Section 4.2.3, the 
connection is by no means straightforward, since, in the present scheme, decoherence 
affects the electronic levels only indirectly. 
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Fig. 17. Maximally achievable concurrence c^ ax , as a function of the number of ions, 
for reservoir coupling T/v = (circles), 1 x 10 -4 (squares), 2 x 10~ 4 (diamonds), 
3 x 10 -4 (triangles) and 4 x 10 -4 (stars). The bold line shows the exact growth of 
the multipartite concurrence as a function of N for GHZ-like states, which satu- 
rates at \/2. This growth is not monotonic anymore when motional decoherence is 
considered, indicating that environment effects become more and more detrimental 
for growing system size. 

0.5 i 1 1 1 1 




Fig. 18. Difference A c = cjy — c^ ax between the ideal and the maximally achievable 
concurrence, as a function of the system size N. Squares, diamonds, triangles and 
stars represent reservoir coupling strengths Y jv = 1 x 10~ 4 , 2 x 10~ 4 , 3 x 10 and 
4 x 10~ 4 , respectively. Consistently with Fig. 17, the entanglement loss during the 
preparation process increases with N. 
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bounds for the various approximations? Are there classes of entangled states 
which cannot be detected by our lower bounds? Is there a suitable generaliza- 
tion of our characterisation of multipartite entanglement through projectors 
on antisymmetric subspaces such as to have a complete description of mul- 
tipartite quantum correlations? And is there a generalisation for continuous 
variable systems? 

And on the physical side: Given the now possible monitoring of the time evolu- 
tion of entanglement under environment coupling, how are the entanglement 
decay rates encoded in the coherences of the original density matrix? How 
do these decay rates depend on i) the initial state and ii) the environment 
coupling operators and strengths? And, perhaps most importantly: Since the 
Hilbert size dimension increases exponentially with increasing systems size 
(i.e. increasing number of system components), quantum state tomography 
will become inoperational to quantify the state's entanglement properties. 
Which experimental observables are then best suited to fulfill this task? Or, 
in other words, which are the robust dynamical observables which exhibit a 
clear experimental signature of the various types of multipartite entanglement? 

Hence, despite the considerable progresses on which we have reported here, 
a panoply of challenging open questions awaits solution, what simply reflects 
the potential - on the fundamental level as well as on the level of applications 
(from decoherence control to quantum computation) - of entanglement theory, 
a field still in its infancy, at the emerging interface of mathematical physics 
and experiments. 
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